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Abstract 

We study the integrable structure of the Dirichlet boundary problem in two 
dimensions and extend the approach to the case of planar multiply-connected do- 
mains. The solution to the Dirichlet boundary problem in multiply-connected case 
is given through a quasiclassical tau-function, which generalizes the tau-function 
of the dispersionless Toda hierarchy. It is shown to obey an infinite hierarchy of 
Hirota-like equations which directly follow from properties of the Dirichlet Green 
function and from the Fay identities. The relation to multi-support solutions of 
matrix models is briefly discussed. 
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1 Introduction 



The Dirichlet boundary problem pQ is to reconstruct a harmonic function in a bounded 
domain from its values on the boundary. Remarkably, this standard problem of complex 
analysis, related however to string theory and matrix models, possesses a hidden inte- 
grable structure [2] , which we clarify further in this paper. It turns out that variation of a 
solution to the Dirichlet problem under variation of the domain is described by an infinite 
hierarchy of non-linear partial differential equations known (in the simply-connected case) 
as dispersionless Toda hierarchy. It is a particular example of the universal hierarchy of 
Whitham equations introduced in |3J E] . 

The quasiclassical tau-function or, more precisely, its logarithm F, is the main new 
object associated with a family of domains in the plane. Any domain in the complex plane 
with sufficiently smooth boundary can be parameterized by its moments with respect to 
a basis of harmonic functions. The F-function is a function of the full infinite set of 
the moments. The first order derivatives of F are then moments of the complementary 
domain. This gives a formal solution to the inverse potential problem, considered for 
the simply-connected case in 10] - The second order derivatives are coefficients of the 
Taylor expansion of the Dirichlet Green function and therefore they solve the Dirichlet 
boundary problem. These coefficients are constrained by infinite number of universal (i.e. 
domain-independent) relations which, unified in a generating form, just constitute the 
dispersionless Hirota equations. For the third order derivatives (their role in problems 
of complex analysis is not yet quite clear) there is a nice "residue formula" which allows 
one to prove (7j that F obeys the WDVV equations. 

Below we are going to demonstrate that for planar multiply-connected domains the 
solution to the Dirichlet boundary problem can be performed in a similar way. Specifi- 
cally, we consider domains which are obtained by cutting several "holes" in the complex 
plane. Boundaries of the holes are assumed to be smooth simple non-intersecting curves. 
In this case, the complete set of independent variables can be again identified with the 
set of harmonic moments. However, a choice of the proper basis of harmonic functions 
in a multiply-connected domain becomes crucial for our approach. It turns out that the 
Laurent polynomials which were used in the simply-connected case should be replaced by 
the basis analogous to the one introduced in jH] - a "global" generalization of the Laurent 
basis for algebraic curves of arbitrary genus. The basis has to be also enlarged to include 
harmonic functions with multi- valued analytic part. This results in an additional finite 
set of extra variables. We construct the F-function and prove that its second derivatives 
satisfy non-linear relations, which generalize the Hirota equations of the dispersionless 
Toda hierarchy. These relations are derived from the Fay identities 9J for the Riemann 
theta functions on the Jacobian of Riemann surface obtained as the Schottky double of 
the plane with given holes. 

We note that extra variables, specific for the multiply-connected case, can be chosen 
in different ways and possess different geometric interpretations, depending on the choice 
of basis of homologically non-trivial cycles on the Schottky double. The corresponding 
F-functions are shown to be connected by a duality transformation - a (partial) Legendre 
transform, with the generalized Hirota relations being the same. 

Now let us give a bit more expanded description of the Dirichlet problem in planar 
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domains. Let D c be a domain in the complex plane bounded by one or several non- 
intersecting smooth curves. It will be convenient to realize D c as a complement to another 
domain D, having one or more connected components, and to consider the Dirichlet 
problem in D c : to find a harmonic function u(z) in D c such that it is continuous up to 
the boundary, <9D C , and equals a given function «o(0 on the boundary. The problem has 
a unique solution written in terms of the Dirichlet Green function G(z,£): 

«(*) = -^-/ u (0d n G(z,0m (1.1) 

2,71 JdD c 

where d n is the normal derivative on the boundary with respect to the second variable, 
the normal vector n is directed inward D c , and \d£\ := dl(£) is an infinitesimal element 
of the length of the boundary <9D C . 

The main object to study is, therefore, the Dirichlet Green function. It is uniquely 
determined by the following properties pQ: 

(Gl) The function G(z, z') is symmetric and harmonic everywhere in D c (including oo 
if D c 3 oo) in both arguments except z = z' where G(z, z') = log \z — z'\ + . . . as 
z — > z'\ 

(G2) G(z, z') = if any one of the variables z, z' belongs to the boundary <9D C . 

Note that the definition implies that G(z,z') < inside D c . In particular, d n G(z,£) is 
strictly negative for all £ G <9D C . 

If D c is simply-connected (note that we assume oo G D c ), i.e., the boundary has only 
one component, the Dirichlet problem is equivalent to finding a bijective conformal map 
from D c onto the complement to unit disk or any other reference domain for which the 
Green function is known explicitly. Such bijective conformal map w(z) exists by virtue 
of the Riemann mapping theorem, then 

w(z) — w(z') 



G(z,z') = log 



1.2 



w(z)w(z') — 1 

where bar means complex conjugation. It connects the Green function at two points with 
the conformal map normalized at some third point (say at z = oo: w(oo) = oo). It is 
this formula which allows one to derive the Hirota equations for the tau-function of the 
Dirichlet problem in the most economic and transparent way [2] (see also sect. El below). 

For multiply-connected domains, formulas of this type based on conformal maps do 
not really exist. In general, there is no canonical choice of the reference domain, moreover, 
the shape of a reference domain depends on D c itself. In fact, as we demonstrate in 
the paper, the correct extension of (jl.2|) needed for derivation of the generalized Hirota 
equations follows from a different direction which is no longer explicitly related to bijective 
conformal maps. Namely, logarithm of the conformal map log w(z) should be replaced 
now by the Abel map from the Schottky double of D c to the Jacobi variety of this 
Riemann surface, and the rational function under the logarithm in (jl.2j) is substituted 
by ratio of the prime forms or Riemann theta-functions. 

We show that the Green function of multiply-connected domains admits a represen- 
tation through the logarithm of the tau-function of the form 

- - ± + \v{z)V{z')F (1.3) 



G(*,*')=log 
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Figure 1: A "pictorial" derivation of the Hadamard formula. We consider a small de- 
formation of the domain, with the new boundary being depicted by the dashed line. 
According to (G2) the Dirichlet Green function vanishes, G(z,£) = 0, if £ belongs to 
the old boundary. Then the variation 6G(z,£) is simply equal to the new value, i.e. in 
the leading order 5G(z,£) = —5n(£)d n G(z,£). Now notice that 8G(z,£) is a harmonic 
function (the logarithmic singularity cancels since it is the same for both old and new 
functions) with the boundary value —5n(£)d n G(z, £). Applying M.l]) one obtains A1.4}) . 
The argument is the same for both simply- connected and multiply-connected domains. 

Here V(z) is certain vector field on the moduli space of boundary curves, therefore it 
can be represented as a (first-order) differential operator w.r.t. harmonic moments with 
constant (in moduli) coefficients depending, however, on the point 2 as a parameter. 

In this paper we also obtain similar formulas for the harmonic measures of the bound- 
ary components and for the Abel map. A combination of these formulas with the Fay 
identities yields the generalized Hirota-like equations for the tau-function F. 

Our main tool is the Hadamard variational formula ^U] which gives variation of the 
Dirichlet Green function under small deformations of the domain in terms of the Green 
function itself: 

5G(z,z') = ±i d n G(z,Z)d n G(z',Z)6n(0\dt\. (1.4) 

2tt JdD c 

Here 5n(£) is the normal displacement (with sign!) of the boundary under the deforma- 
tion, counted along the normal vector at the boundary point £. It was shown in [2] that 
this remarkable formula is a key to all integrable properties of the Dirichlet problem. An 
extremely simple "pictorial" derivation of the formula (|1.4jl is presented in fig. ^ 

We start with a brief recollection of the results for the simply-connected case in sect. |2] 
However, instead of "bump" deformations used in j2] we work here with their rigorously 
defined versions - a family of infinitesimal deformations which we call elementary ones. 
This approach is basically motivated by the theory of interface dynamics in viscous fluids, 
which is known to be closely connected with the formalism developed in [2] and in the 
present paper (see jHj for details). 

In sect. E] we introduce local coordinates in the space of planar multiply-connected do- 
mains and express the elementary deformations in these coordinates. Using the Hadamard 
formula, we then observe remarkable symmetry or "zero-curvature" relations which con- 
nect elementary deformations of the Green function and harmonic measures. The ex- 
istence of the tau-function and the formula (jl.3j) for the Green function directly follow 
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from these relations. In sect. 0]we make a Legendre transform to another set of local 
coordinates in the space of algebraic multiply-connected domains, which is in a sense 
dual to the original one. In these coordinates, eq. (jl.3|) gives another version of the 
Green function which solves the so-called modified Dirichlet problem. We also discuss 
the relation to multi-support solutions of matrix models in the planar large N limit. 

In sect. 5 we combine the results outlined above with the representation of the Green 
function in terms of the prime form on the Schottky double. This allows us to obtain an 
infinite system of partial differential equations on the tau-function which generalize the 
dispersionless Hirota equations. 



2 The Dirichlet problem for simply-connected do- 
mains and dispersionless Hirota equations 

In this section we rederive the results from [2] for the simply-connected case in a slightly 
different manner, more suitable for further generalizations. At the same time we show 
that the results of |2j obtained for analytic curves can be easily extended to the smooth 
case. 

Let D be a connected domain in the complex plane bounded by a simple smooth curve. 
We consider the exterior Dirichlet problem in D c = C \ D which is the complement of 
D in the whole (extended) complex plane. Without loss of generality, we assume that D 
is compact and contains the point z = 0. Then D c is a simply-connected domain on the 
Riemann sphere containing oo. 



Harmonic moments and deformations of the boundary. Let t k be moments of 
the domain D c = C \ D defined with respect to the harmonic functions {z~ k /k}, k > 0: 

t k = -\ f z- k d 2 z, k = l,2,... (2.1) 

TCK JD C 

and {i k } be the complex conjugate moments, i.e. i k = — ^ J Dc d 2 zz~ k . The Stokes 
formula represents the harmonic moments as contour integrals 

** = 7tW ^ k zdz ( 2 - 2 ) 

llXlk JdD 

providing, in particular, a regularization of possibly divergent integrals (|2.1jl . Besides, 
we denote by to the area (divided by n) of the domain D: 

t Q = - [ d 2 z (2.3) 

7T JD 

The harmonic moments of D c are coefficients of the Taylor expansion of the potential 

$M = -- f \og\z- z'\d 2 z' (2.4) 

7T JD 

induced by the domain D filled by two-dimensional Coulomb charges with the uniform 
density p = —1. Clearly, d z dg&(z) = — 1 if z G D and vanishes otherwise, so around the 
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origin (recall that D 3 0) the potential equals to — \z\ 2 plus a harmonic function, i.e. 

S(z) - $(0) = -\z\ 2 + £ + (2 5) 

fe>i 

and one can verify that t& are just given by (|2.1jl . 

For analytic boundary curves, one may introduce the Schwarz function associated 
with the curve. The function 

1 r d 2 z' 



71 JD Z 

is continuous across the boundary and holomorphic for zeD c while for z G D the function 
d z $ + z is holomorphic. If the boundary is an analytic curve, both these functions can 
be analytically continued outside the regions where they were originally defined, and, 
therefore, there exists a function, S(z), analytic in some strip- like neighborhood of the 
boundary contour, such that S(z) = z on the contour. In other words, S(z) is the analytic 
continuation of z away from the boundary contour, this function completely determines 
the shape of the boundary and is called the Schwarz function jllj . In general we are 
going to work with smooth curves, not necessarily analytic, when the Schwarz function 
does not exist as an analytic function. Nevertheless, it appears to be useful below to 
define the class of boundary contours with nice algebro-geometric properties. 

The basic fact of the theory of deformations of closed smooth curves is that the (in 
general complex) moments {tk,ik} = {t±k} supplemented by the real variable to form a 
set of local coordinates in the "moduli space" of smooth closed curves [T2] (see also [T5|). 

Important remark. This means that: (a) under any small deformation of the domain the 
set t = {to,t±k} is subject to a small change; (b) on the space of smooth closed curves 
there exist vector fields d tk such that d tk t n = 5k n , which are represented in terms of 
infinitesimal normal displacements of the boundary that change = Re t^ or = Im t k 
keeping all the other moments fixed; (c) the corresponding infinitesimal displacements 
can be locally integrated. The latter means that for each domain D c with moments 
{to,t±k} and for an arbitrary integer N there exist constants e m , \m\ < N, such that for 
any set {t' ,t' ±k } with \t' m — t m \ < e m , m < N, t' m = t m , \m\ > N, in the neighborhood 
of D c there is a unique domain with the moments {t' Q ,t± k }. We adopt this restricted 
notion of the local coordinates throughout the paper. It would be very interesting to find 
conditions on the infinite sets for the corresponding rectangles to form an open set 
in an infinite-dimensional variety of smooth curves. We plane to address this problem 
elsewhere. 

Let us present a proof of this statement which later will be easily adjusted to the 
case of multiply-connected domains. At the same time this proof allows one to derive 
a deformation of the domain with respect to the variables Suppose there is a one- 
parametric deformation D(t) (with some real parameter t) of D = D(0) such that all 
are preserved: d t tk = 0, k > 0. Let us prove that such a deformation is trivial. The proof 
is based on two key observations: 

• The difference of the boundary values d t C ± (()d( of the derivative of the Cauchy 
integral 

C(z)dz = ^i M- (2.6) 
2m J 3D Q — z 
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is purely imaginary differential on the boundary of D. 

Indeed, let ({a,t) be a parameterization of the curve dD(t). Denote the value of 
the differential (|2.fi|) by C~{z)dz for zgD c and by C + (z)dz for z G D. Taking the 
t-derivative of (|2.fij) and integrating by parts one gets 



2,711 JdD y C — z [C — z) 1 ) 



da 



dz r ( ( t ( a - CaCt 
2m JdD \ C — z 

Hence, 

(d t C+(C) - d t C-(()) d( = d t (d( - d t (d( = Mm (d t (d( 
is indeed purely imaginary. 



(2.7) 



• If a t- deformation preserves all the moments t^, k > 0, i/ie differential 
d t (d( — d t (d( extends to a holomorphic differential in D c . 

If \z\ < |C| for all C e 9D, then we can expand: 

f) / d? 00 r - \ 00 

9 ' c+ ^ dz = mWiE zt L^ d( ) =EM^)^ = o (2.8) 

and, since C + is analytic in D, we conclude that <9iC + = 0. The expression d t (d( — 
d t (d( is the boundary value of the differential —d t C~(z)dz which has at most simple 
pole at the infinity and holomorphic everywhere else in D c . The equality 

dth = ^-i (d t (d(-d t (d() = 
2m JdD 

then implies that the residue at z = oo vanishes, therefore d t C~{z)dz is holomor- 
phic. 

Any holomorphic differential which is purely imaginary along the boundary of a 
simply-connected domain must be zero in this domain. Indeed, the real part of the 
harmonic continuation of the integral of this differential is a harmonic function with a 
constant boundary value. Such a function must be constant by virtue of the uniqueness 
of the solution to the Dirichlet problem. Another proof relies on the Schwarz symmetry 
principle and the standard Schottky double construction (see the next section for de- 
tails). Consider the compact Riemann surface obtained by attaching to D c its complex 
conjugated copy along the boundary. Since d t C~dz is imaginary along the boundary, 
we conclude, from the Schwarz symmetry principle, that dtC~dz extends to a globally 
defined holomorphic differential on this compact Riemann surface, which has genus zero. 
Therefore, such a differential is equal to zero. Hence we conclude that d t (d( — d t (d( = 0. 
This means that the vector d t ( is tangent to the boundary. Without loss of generality 
we can always assume that a parameterization of dD(t) is chosen so that d t C,{a,i) is nor- 
mal to the boundary. Thus, the t-deformation of the boundary preserving all harmonic 
moments is trivial. 

The fact that the set of harmonic moments is not overcomplete follows from the 
explicit construction of vector fields in the space of domains that changes any harmonic 
moment keeping all the others fixed (see below). 
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Figure 2: The elementary deformation with the base point z. 

Elementary deformations and the operator V(z). Fix a point z G D c and consider 
a special infinitesimal deformation of the domain such that the normal displacement of the 
boundary is proportional to the gradient of the Green function G(z,£) at the boundary- 
point (fig. El): 

6n{Z) = ~d n G{z,£) (2.9) 

For any sufficiently smooth initial boundary this deformation is well-defined as e — > 0. We 
call infinitesimal deformations from this family, parametrized by 2 6 D c , the elementary 
deformations. The point z is refered to as the base point of the deformation. Note 
that since d n G < (see the remark after the definition of the Green function in the 
Introduction), 5n for the elementary deformations is either strictly positive or strictly 
negative depending of the sign of e. 

Let 5 Z be variation of any quantity under the elementary deformation with the base 
point z. It is easy to see that 5 z t = e, 5 z tk = ez~ k /k. Indeed, 

S z t = - fSn(£M\ = ~ id n G(z,0m = e 

71 J Lix J 

(2.10) 

6ut k = ±fc"wtM\ = -^/r*a.G(*,oidei = 

by virtue of the Dirichlet formula (|1.1J1 . Note that the elementary deformation with the 
base point at oo keeps all moments except to fixed. Therefore, the deformation which 
changes only to is given by 5n(£) = — f<9 n G(oo, £). 

Now we can explicitly define the deformations that change only either Xk = Ret^ or 
Uk = Imtfe keeping all other moments fixed. As is clear from (J2.1(J)) . the corresponding 
8n(£) is given by the real or imaginary part of normal derivative of the function 

Hk(0 = ^-<f z k d z G{z,i)dz (2.11) 

L1XI Joo 

at the boundary. Here the contour integral goes around infinity. Namely, the normal 
displacements 5n(£) = eRe (d n Hk(£)) and 6n(£) = elm (d n Hk(£)) change the real and 
imaginary part of tk by ±e respectively keeping all other moments fixed. 
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These deformations allow one to introduce the vector fields 

d d d 
dt ' dx k ' dy k 

in the space of domains which are locally well-defined. Existence of such vector fields 
means that the variables are independent. For k > it is more convenient to use their 
linear combinations 

— 1 ( 9 i 9 \ — 1 ( 9 | j 9 \ 

dt k 2\dx k dyk)' di k 2 \dx k dy k ) 

which span the complexified tangent space to the space of simply-connected domains 
(with fixed area to)- If X is any functional of our domain locally representable as a 
function of harmonic moments, X = X(t), the vector fields dt , dt k , di k can be understood 
as partial derivatives acting to the function X(t). 

Consider the variation 5 Z X of a functional X — X(t) under the elementary deforma- 
tion with the base point z. In the leading order in e we have: 

dX 

<^ X = E ^7" 5 ^k = tV{z)X (2.12) 
k atk 

where the differential operator V(z) is given by 

V(*) = d t0 + + (2.13) 

The right hand side suggests that for functionals X such that the series V(z)X converges 
everywhere in D c up to the boundary 5 Z X is a harmonic function of the base point z. 

Note that in j2] we have used the "bump" deformation and continued it harmonically 
into D c . In fact, it was the elementary deformation (JTHIJ) 5 Z oc § \d^\d n G{z, £,)5 hump (£) 
that was really used. The "bump" deformation should be understood as a (carefully 
taken) limit of 5 Z when the point z tends to the boundary <9D C . 



The Hadamard formula as integrability condition. Variation of the Green func- 
tion under small deformations of the domain is known due to Hadamard, see eq. (|1.4|) . To 
find how the Green function changes under small variations of the harmonic moments, 
we fix three points a,b,c G C \ D and compute 5 c G(a, b) by means of the Hadamard 
formula ([1.4)1 . Using (|2.12j) . one can identify the result with the action of the vector field 
V(c) on the Green function: 

V(c)G(a,6) = --^/ d n G(a,0d n G(b,0d n G(c,0m (2.14) 

4-7T JdD 

Remarkably, the r.h.s. of (|2.14|) is symmetric in all three arguments, i.e. 

V(a)G(6, c) = V(6)G(c, a) = V(c)G(a, b) (2.15) 

This is the key relation which allows one to represent the Dirichlet problem as an in- 
tegrable hierarchy of non- linear differential equations 2 , (12.15(1 being the integrability 
condition of the hierarchy. 
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It follows from (j2.15j) (see [2] for details) that there exists a function F = F(t) such 



that 



G(z,z')= log 



+ -V(z)V(z')F 



(2.16) 



We note that existence of such a representation of the Green function was first conjectured 
by Takhtajan. For the simply-connected case, this formula was obtained in (see also 
[T3] for a detailed proof and discussion). The function F is (logarithm of) the tau-function 
of the integrable hierarchy. In ^3] it was called the tau-function of the (real analytic) 
curves - the boundary contours <9D or <9D C . 



Dispersionless Hirota equations. Combining f)2.16j) and (jl.2j) . we obtain the relation 



log 



w(z) — w(z') 



w(z)w(z') — 1 



2 


1 1 


= log 






z z' 



+ V(z)V(z')F 



(2.17) 



which implies an infinite hierarchy of differential equations on the function F. It is 
convenient to normalize the conformal map w(z) by the conditions that w(oo) = oo and 
d z w(oc) is real, so that 

w(z) = - + OCl) asz^oo (2.18) 
r 

where the real number r = lim^oo dz/dw(z) is called the (external) conformal radius 
of the domain D (equivalently, it can be defined through the Green function as logr = 
]im. z ^ OQ (G(z, oo) + log \z\), see [15J). Then, tending z' — > oo in (|2.17)L one gets 

log \w(z)\ 2 = log \z\ 2 - d t0 V(z)F (2.19) 

The limit z — > oo of this equality yields a simple formula for the conformal radius: 

logr 2 = «9 2 F (2.20) 

Let us now separate holomorphic and antiholomorphic parts of these equations, in- 
troducing the holomorphic and antiholomorphic parts of the operator V(z) (J2.13I) : 



k>l 



fc>i K 



(2.21) 



Rewrite (|2.17J) in the form 



( W {z)-, { zi)\_ n__i 

\w{z)w{z>) -l) B \z z> 



d t0 + D(z) ) V(z')F 



-log 



w(z) — w(z') 
w{z')w(z) 



The l.h.s. is a holomorphic function of z while the r.h.s. is antiholomorphic. Therefore, 
both are equal to a z- independent term which can be found from the limit z — > oo. As 
a result, we obtain the equation 



log 



w(z) — w(z') 
jw(z) - (w(z'))-\ 



log 1-- +D{z)V{z')F 



(2.22) 
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which, as z' — > oo, turns into the formula for the conformal map w(z): 



\ogw(z) = log* - l -d 2 tQ F - d t0 D(z)F (2.23) 



(here we also used (|2.20Jl ). Proceeding in a similar way, one can rearrange (J2.22)) in order 
to write it separately for holomorphic and antiholomorphic parts in z'\ 

[og m W -.,(,') = _^ F + D{Z)D{Z , )F (2 . 24) 

- log ( 1 - , 1 ) = D(z)D(z')F (2.25) 
\ w{z)w(z') j 

Writing down eqs. (|2.24|) for the pairs of points (a, b) , (6, c) and (c, a) and summing 
up the exponentials of the both sides of each equation one arrives at the relation 

( fl _ b ) e D(a)D(b)F + (6 _ c)e D(b)D { c)F + ^ _ a)e D(c)D(a)F = g 

which is the dispersionless Hirota equation (for the KP part of the two-dimensional Toda 
lattice hierarchy) written in the symmetric form. This equation can be regarded as a 
very degenerate case of the trisecant Fay identity 0. It encodes the algebraic relations 
between the second order derivatives of the function F. As c — > oo, we get these relations 
in a more explicit but less symmetric form: 

1 _ e D(a)D {b )F = D(a) - D(b) d F (2>27) 
a — b 

which makes it clear that the totality of second derivatives := d^d^F are expressed 
through the derivatives with one of the indices put equal to unity. 

More general equations of the dispersionless Toda hierarchy obtained in a similar way 
by combining eqs. (J2.2Hj) . ()2.24|) and (|2.25j) include derivatives w.r.t. to and tk- 

( a _ b)e D(a)D(b)F = ae -9 t0 D(a)F _ ^-d^F ^ 

1 _ e -D(z)D(z)F = ]_ e dt V(z)F (2.29) 

ZZ 

These equations allow one to express the second derivatives d tm d tn F, d tm di n F with 
m, n > 1 through the derivatives d to d tk F, d to di k F. In particular, the dispersionless 
Toda equation, 

d tl d h F = e d ?o F (2.30) 
which follows from 1)2.29)1 as z — * oo, expresses d^d^F through d 2 Q F. 

For a comprehensive exposition of Hirota equations for dispersionless KP and Toda 
hierarchies we refer the reader to [To*) IT7] . 



11 



Integral representation of the tau-function. Eq. (j2.16|) allows one to obtain a rep- 
resentation of the tau-function as a double integral over the domain D. Set &(z) : = 
V(z)F. One is able to determine this function via its variation under the elementary 
deformation: 

S a $( z ) = -2elog|a" 1 - z^\ + 2eG(a,z) (2.31) 

which is read from eq. (j2.1tij) by virtue of ()2.12|) . This allows one to identify $ with the 
"modified potential" $(2;) = $(z) - $(0) + t log \z\ 2 , where $ is given by ()2.4j) . Thus we 
can write 



V(z)F = $(z) = -- I log|^ 1 -rV 2 C = ^o + 2Re^^- fc (2.32) 

7t JD k>Q fc 



The last equality is to be understood as the Taylor expansion around infinity. The 
coefficients Vk are moments of the interior domain (the "dual" harmonic moments) defined 
as 

v k = - [ z k d 2 z (k > 0) , v = -$(0) = - / log \z\d 2 z (2.33) 

It JD It JD 

From ()2.32|) it is clear that 

v k = d tk F, k>0 (2.34) 

i.e., the moments of the complementary domain D (the "dual" moments) are completely 
determined by the function F of harmonic moments of D c . 

In a similar manner, one arrives at the integral representation of the tau-function. 
Comparing ()2.32j) with ()2.31|) one can easily notice that the meaning of the elementary 
deformation 5^ or the operator V(£) formally applied at the boundary point £ G 3D 
(where G(z,C,) = 0) is attaching a "small piece" to the integral over the domain D (the 
"bump" operator from [2 ). Using this fact and interpreting ()2.32|) as a variation 5 Z F we 
arrive at the double-integral representation of the tau-function 

F = -— I I log \z~ x - CV 2 ^ 2 C (2-35) 

or 



7T 2 JDJD 



F = — [ $(z)d 2 z=— [ ($(z) - 2$(0)) d 2 z (2.36) 

lit JD lit JD 

As we see below, the main formulas from this paragraph remain intact in the multiply- 
connected case. 



3 The Dirichlet problem and the tau-function in the 
multiply-connected case 

Let now D a , a — 0, 1, . . . , g, be a collection of g + 1 non-intersecting bounded connected 
domains in the complex plane with smooth boundaries dD a . Set D = U^qDq,, so that 
the complement D c = C \ D becomes a multiply-connected unbounded domain in the 
complex plane (see fig. EJ). Let be the boundary curves. They are assumed to be 
positively oriented as boundaries of D c , so that U 9 a=0 b a = dD c while b a = — dD a has the 
clockwise orientation. 
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Figure 3: A multiply-connected domain D c = C \ D for g = 3. The domain D = U«=o D a 
consists of g + 1 = 4 disconnected parts D a with the boundaries b a . To define the 
complete set of harmonic moments, we also need the auxiliary points z a G D a which 
should be always located inside the corresponding domains. 

Comparing to the simply-connected case, nothing is changed in posing the standard 
Dirichlet problem. The definition of the Green function and the formula (jl.lj) for the 
solution of the Dirichlet problem through the Green function remain to be the same. 

A difference is in the nature of harmonic functions. Any harmonic function is still 
the real part of an analytic function but in the multiply-connected case these analytic 
funstions are not necessarily single- valued (only their real parts have to be single- valued) . 
In other words, the harmonic functions may have non-zero "periods" over non-trivial 
cycles 1 . In our case, the non-trivial cycles are the boundary curves b Q . In general, the 
Green function has non-zero "periods" over all boundary contours. Hence it is natural 
to introduce new objects, specific to the multiply-connected case, which are defined as 
"periods" of the Green function. 

First, the harmonic measure w a (z) of the boundary component b a is the harmonic 
function in D c such that it is equal to unity on b a and vanishes on the other boundary 
curves. Thus the harmonic measure is the solution to the particular Dirichlet problem. 
From the general formula (jl.lj) we conclude that 

u a (z) = -—<f d n G(z,()\d(\, a = l,...,g (3.1) 

Z7T Jb a 

so the harmonic measure is the period of the Green function w.r.t. one of its arguments. 
From the maximum principle for harmonic functions it follows that < u a (z) < 1 in 
internal points. Obviously, J2 9 a =o u a (z) = 1. In what follows we consider the linear 
independent functions u a (z) with a = 1, . . . , g. 

1 Here and below by "periods" of a harmonic function / we mean the integrals § d n f dl over non-trivial 
cycles. 
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Further, taking "periods" in the remaining variable, we define 



i d n u a (()\d(\, a,p = l,...,g 



(3.2) 



The matrix Q a p is known to be symmetric, non-degenerate and positively-definite. It will 
be clear below that the matrix T a p = mVL a p can be identified with the matrix of periods 
of holomorphic differentials on the Schottky double of the domain D c (see formula (j3.5|0 . 
For brevity, we refer to both T a p and VL a p as period matrices. 

For the harmonic measure and the period matrix there are variational formulas similar 
to the Hadamard formula (|1.4|) . They can be derived either by a direct variation of (|3.1|) 
and ()3.2j) using the Hadamard formula or (much easier) by a "pictorial" argument like 
in fig. Q The formulas are: 



The Schottky double. It is customary to associate with a planar multiply-connected 
domain its Schottky double (see, e.g., ^H], Ch. 2.2), a compact Riemann surface without 
boundary endowed with antiholomorpic involution, the boundary of the initial domain 
being the set of the fixed points of the involution. The Schottky double of the multiply- 
connected domain D c can be thought of as two copies of D c ( "upper" and "lower" sheets 
of the double) glued along the boundaries U 9 a=0 b a = <9D C , with points at infinity added 
(oo and oo). In this set-up the holomorphic coordinate on the upper sheet is z inherited 
from D c , while the holomorphic coordinate 2 on the other sheet is z. The Schottky double 
of D c with two infinities added is a compact Riemann surface £ of genus g = #{D Q } — 1. 
A meromorphic function on the double is a pair of meromorphic functions /, / on D c such 
that f(z) = f(z) on the boundary. Similarly, a meromorphic differential on the double is 
a pair of meromorphic differentials f(z)dz and f(z)dz such that f(z)dz = f(z)dz along 
the boundary curves. 

On the double, one may choose a canonical basis of cycles. We fix the b-cycles to be 
just the boundaries of the holes b a for a = 1, . . . , g. Note that regarded as the oriented 
boundaries of D c (not D) they have the clockwise orientation. The a a -cycle connects the 
a-th hole with the 0-th one. To be more precise, fix points £ a on the boundaries, then 
the a a cycle starts from £ , goes to £ Q on the "upper" (holomorphic) sheet of the double 
and goes back the same way on the "lower" sheet, where the holomorphic coordinate is 
z, see fig. EJ 

Being harmonic, uo a can be represented as the real part of a holomorphic function: 



where W a (z) are holomorphic multivalued functions in D c . The differentials dW a are 
holomorphic in D c and purely imaginary on all boundary contours. So they can be 

2 More precisely, the proper coordinates should be l/z (and l/z), which have first order zeros instead 
of poles at z = oo (and z = do). 




(3.3) 




(3.4) 



OJa(z) = W a {z) + W a ( Z ) 
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Figure 4: The domain D c with the a a -cycle, going one way along the "upper sheet" and 
back along the "lower sheet" of the Schottky double of D c . For such choice one clearly 
gets the intersection form a a o bp = 5 a /3 for a, (3 = 1, . . . , g. 



extended holomorphically to the lower sheet as —dW a (z). In fact this is the canonically 
normalized basis of holomorphic differentials on the double: according to the definitions, 

/ dW p = t dWp(z) + / € ° (-dW^zf) = 2Re [** dWp(z) = Wf>(£ a )-M£o) = <W 

Then the matrix of b-periods of these differentials reads 

T Q/3 = f dWp = - % - i d n updl = mVt al3 (3.5) 

i.e. the period matrix T a @ of the Schottky double S is a purely imaginary non-degenerate 
matrix with positively definite imaginary part ftVt a p (|3.2j) . 



Harmonic moments of multiply-connected domains. One may still use harmonic 
moments to characterize the shape of a multiply-connected domain. However, the set 
of harmonic functions should be extended by adding functions with poles in any hole 
(not only in D as before) and functions whose holomorphic parts are not single-valued. 
To specify them, let us mark points z a e D a , one in each hole (see fig. EI). Without 
loss of generality, it is convenient to put zq = 0. Then one may consider single- valued 



\z - z n 



and harmonic functions log 



with 



analytic functions in D c of the form 
multi- valued analytic part. 

The arguments almost identical to the ones used in the simply-connected case show 



that the parameters t , M n 
M„ a ■- 



7T JD' 



where as in ()2.3|) t = Area(D)/7r, 
(z — z a )~ n d 2 z, a = 0, 1, . . . , g, n>l 



together with their complex conjugate, and 

1 



7T 



/ log 






Z 



d z, 



a 



(3.6) 



(3.7) 
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uniquely define D c , i.e. any deformation preserving these parameters is trivial. Note that 
the extra moments <j) a are essentially the values of the potential (|2.4|) at the points z a 

<f) a = $(0) -Q(z a ) - \z a \ 2 (3.8) 

A crucial step for what follows is the change of variables from M n>a to the variables 
Tfc which are finite linear combinations of the M n ^s. They can be directly defined as 
moments with respect to new basis of functions: 

= to, n = 7T— / A k (z)zdz = -- / d 2 zA k (z), k > (3.9) 
2m JdD it Jd c 

The functions A k (z) are analogous to the Laurent-Fourier type basis on Riemann surfaces 
introduced in jH] . They are explicitly defined by the following formulas (here the indices 
a and (3 are understood modulo g + 1): 

a-l 

A m{g+l)+a = R~ m (z) H(z- zp)- 1 

0=0 

R(z) = f[(z - Zp ) 

(3=0 

In a neighbourhood of infinity A k (z) = z~ k + 0(z^ k ^ 1 ). Any analytic function in D c 
vanishing at infinity can be represented as a linear combination of A k which is convergent 
in domains such that |i2(z)| > const. In the case of one hole (g = 0) the formulas (|3.10|) 
give the basis used in the previous section: A k = z~ k . Note that Aq — 1, A\ — 1/z, 
therefore r = t and t\ = M 1)0 = t\. 



(3.10) 



Local coordinates in the space of multiply-connected domains. Now we are 

going to prove that the parameters Tfc, 4> a can be treated as local coordinates in the 
space of multiply-connected domains. (Here we use the same restricted notion of the 
local coordinates, as in the simply connected case (see Remark in the Section 2)). 

It is instructive and simpler first to prove this for another choice of parameters. 
Instead of <j) a one may use the areas of the holes 

s a = Area = - / d 2z = ^—-i > zdz, a = l,...,g (3.11) 

7T 7T JD a 2m JdD a 

In order to prove that any deformation that preserves r k and s a is trivial, we introduce 
the basis of differentials dB k which satisfy the defining "orthonormality" relations 

— / A k dB k > = 5 k ,k' (3-12) 
2m JdD 

for all integer k, k'. It is easy to see that explicitly they are given by: 

dzR m (z) 

dB m(g+1)+a = — — ±-L [[(z- z/3_i) (3.13) 

Z Z 9 (3=0 

where we identify z_i = z g . The existence of a well-defined "dual" basis of differentials 
obeying the orthonormality relation is the key feature of the basis functions A k , which 
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makes r k good local coordinates comparing to the M Ut(X . For the functions (z — z a ) n one 
cannot define the dual basis. 



The summation formulas 

dzd( 



J2d(A n (C)dB n (z), \R(z)\ < \R(0\ 

(3.14) 



dzdC 00 

= - £ d(A_ n (()dB_ n (z), \R(z)\ > \R(()\ 



n=0 



which can be checked directly, allow us to repeat arguments of sect. |2] Indeed, the 
Cauchy integral (|2.6|) . 

C(z)dz = ±f M- (3.15) 

where the integration now goes along all boundary components, defines in each of the 
holes D Q analytic differentials C a (z)dz (analogs of C + (z)dz in the simply-connected case). 
In the complementary domain D c the Cauchy integral still defines the differential C~(z)dz 
holomorphic everywhere in D c except for infinity where it has a simple pole. The differ- 
ence of the boundary values of the Cauchy integral is equal to z: 

C a {z) -C~{z) = z, zedD a 

From equation (J2.7jl . which can be written separately for each contour, it follows that 

• The difference of the boundary values 

(d t C a (0-d t C-(())d( 

of the derivative of the Cauchy integral hS.lfty is, for all a, a purely imaginary 
differential on the boundary b a . 



The expansion ()3.14|) of the Cauchy kernel implies that 

• If a t-deformation preserves all the moments t^, k > 0, then dt(d( — dt(d( extends 
to a holomorphic differential in D c . 

Indeed, since |-R(z)| is small for z close enough to any of the points z a , one can expand 
d t C a (z) for any a as 

d t C a (z)dz = — E dB k d t f A k (()(d() = d t r k dB k (z) (3.16) 
Zm k=i ^ J9D ' k=i 

and conclude that it is identically zero provided d t r k = 0. Hence — d t C~(z)dz is the 
desired extension of d t (d( — d t (d(. It has no pole at infinity due to the equation d t r = 0. 

Using the Schwarz symmetry principle we obtain that dtC~(z)dz extends to a holo- 
morphic differential on the Schottky double. If the variables s a are also preserved under 
the t-deformation, then this holomorphic differential has zero periods along all the cycles 
b a . Therefore, it is identically zero. This completes the proof of the statement that any 
deformation of the domain preserving r k and s a is trivial. 
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In this proof the variables s a were used only at the last moment in order to show 
that the extension of d t C~(z)dz as a holomorphic differential on the Schottky double 
is trivial. The variables <p a can be used in a similar way. Namely, let us show that if 
they are preserved under t-deformation then a a -periods of the extension of d t C~(z)dz 
are trivial, and therefore this extension is identically zero. Indeed, the variable <p a (|3.7|) 
can be represented in the form 

2 /■** r d 2 C 
<P a = —Re / dz / —±- (3.17) 

7T JO JD C Z — Q 

The differential ^ jp c ^ is equal to C a (z)dz for z e D a and (z + C~{z))dz for z G D c . 
Let £o> be the points where the integration path from to z a intersects the boundary 
contours bo, b a . Then 

<j) a = -2 Re ^jf & C°(z)dz + J*" C a (z)dz + \z + C~{z))dz\ (3.18) 

It is shown above that if a t-deformation preserves the variables then all d t C a (z)dz = 0. 
Thus vanishing of the t-derivative d t 4> a = implies 




(3.19) 



The r.h.s of this equation is just the a a -period of the holomorphic extension of the 
differential d t C~(z)dz. 

Let us construct the deformations d% k and d^ k of the boundary that change the real 
or imaginary parts of the variable = Xk + iyki k > 1, keeping all the other moments 
and the variables a fixed. It is convenient to set df k = \{d% k ~ idy k )- The argument is 
similar to the proof of the fact that any deformation that preserves all the variables is 
trivial. 



• Suppose that the deformations d£ and d^ k exist. Then the differential df k (d( — 
df k (d( extends from dD c to the Schottky double S. Its extension is a meromorphic 
differential dflf. with the only pole at the infinity point oo on the upper sheet. In a 
neighborhood of oo it has the form 

dnt(z) = dB k {z) + 0(z~ 2 )dz (3.20) 
The ^-periods of dQf. are equal to 

I dVlf= r dB k (3.21) 

Ja a JO 

First of all, it is clear that the meromorphic differential dQf. on E is uniquely defined by 
its asymptotics at oo and by the normalization (|3.21|) of its a-periods. To deduce these 
properties, we notice, using (j3.16jl . that d Xk C a (z)dz = dB k (z). Therefore, the differential 
d$ k (d( - d£ k C<K extends to D c as 

dn* k = -d Xk C~(z)dz + dB k (3.22) 
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Using the Schwarz symmetry principle we conclude that it extends to the Schottky double 
as a meromorphic differential. Around the two infinities it has the form dQt — dB k + 

0(z~ 2 )dz and dQi _=_ —dB k -\-0(z~ 2 )dz. In the same way one gets that the differential 

_ fc z — >oo 

d^ n C,dC, — d^ n (d( extends to the double as a meromorphic differential dQ^ k , which at the 
two infinities has the form dQt = idB k + 0(z~ 2 )dz and dQt = idB k + 0{z~ 2 )dz 
respectively. Since 

2dQt = dQt k - idQ* k 
the first statement is proven. From d* h <f> a = and (EHTf) . (l3~2~2l it follows that 

= Re f jf & dB k + (-dQt h + dB k ) + jf*" dB k ^j 

Hence 

f dQt k = 2Re / dB k 

Ja a JO 

In the same way one gets 

I dQ^ k = -21m [ Za dB k 

J 3 a JO 

The last two equations are equivalent to ()3.21|) . 

Normal displacement of the boundary that accomplishes the deformations can be 
explicitly found using the following elementary proposition. 

• Let D(t) be a deformation with real parameter t such that the differential 

dQ = d t (d( - d t (d( 

extends to a meromorphic differential dQ globally defined on the Schottky double 
E. Then the corresponding normal displacement of the boundary is proportional to 
normal derivative of Re J z dQ at the boundary point £: 

,Jn(0 = ^5t d n ^Re J* dnj (3.23) 

Conversely, if 6n(£) = ^5td n H(£), where H is a real-valued function such that 
dH = along the boundary contours and d z H is meromorphic in D c then the 
differential d t (d( — d t (d( is meromorphically extendable to the Schottky double as 
2d z Hdz on the upper sheet and —2d z Hdz on the lower sheet. 

In our case normal displacements of the boundary that change x k or y k keeping all the 
other moments and the variables <p a fixed are thus given by 

6n(0 = l -5x k d n (Re f dQ^ , 5n(£) = \by k d n (bb f dfijj (3.24) 

respectively. Note that since the differentials dQ Xk (z), dQ yk (z) (but not dQ k (z)\) are 
purely imaginary on the boundaries, d~ReQ Xk (z) = d~ReQ yk (z) = along each compo- 
nent of the boundary. With formulas (|3.24j) at hand, one can directly verify that these 
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deformations indeed change x k or only and keep fixed all other moments. We leave 
this to the reader. 

In terms of the differential dflf. formulas ()3.24j) acquire the form 

5n(0 = Sx k d n (lie J* dtt{\ , 8n(£) = -5y k d n (im J* dfltj (3.25) 

(cf. (I2.11J) for the simply-connected case). Indeed, taking the real part of 2f4(£) = 
£Xc fc (£) — ^3/ fc (0) we § e t 2Ref2 fc (£) = RelX,. fe (£) + lmfl yk (^). But the normal derivative 
of Im Qy k (£) vanishes since, by virtue of the Cauchy-Riemann identities, it is equal to the 
tangential derivative of the conjugate harmonic function Ref^ fc (£). This proves the first 
formula in ()3.25|) . The second one is proven in a similar way by taking imaginary part 
of 2fi fc (0- 

The construction of the vector fields df o (which changes tq only) and (which changes 
<f) a only) is quite similar and even simpler since the derivative (j3.16|) vanishes. So, we 
present the results without going into details. 



• The deformation df corresponds to the normal displacement 

Sn(C) = --6T d n G(oa,£) 

The differential — (df (d( — df (d() extends from <9D C to the Schottky double E. Its 
extension is a meromorphic third-kind Abelian differential dflo which has simple 
poles at the infinities on the two sheets of the Schottky double (with residues ±1) 
and vanishing ^-periods. 

• The deformation corresponds to the normal displacement 

where u a is the harmonic measure of the boundary component b a (see The 
differential d^QdC, — d^(d( holomorphically extends from <9D C to the Schottky double 
E. Its extension is the canonically normalized holomorphic differential dW a = 
d z uj ol {z)dz on the upper sheet (and dW a = —d g uj a (z)dz on the lower sheet). 



So we see that d% k , d$ k , df and are well-defined vector fields on the space of 
multiply-connected domains. This fact allows us to treat <p a , as local coordinates on 
this space. At this stage it becomes clear why we prefer to use the moments rather 
than Mfc Q . Although the latter are finite linear combinations of the former, they can 
not be treated as local coordinates because the vector fields d/dM kjOC , being in general 
infinite linear combinations of the df , are not well-defined. 



Il-variables. Up to now the roles of the variables s a and <p a have been in some sense 
dual to each other. It is necessary to emphasize that this duality does not go beyond the 
framework of our proof of the statement that the first or the second sets together with 
the variables r k are local coordinates in the space of multiply-connected domains. For 
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analytic boundary curves one can define the Schwarz function, which is a unique function 
analytic in some strip-like neighborhoods of all boundaries such that 

S(z) — z on the boundary curves (3.26) 

Then the variables s a are b-periods of the differential S(z)dz. At the same time, the vari- 
ables 4> a in general can not be identified with periods of this differential (or its extension) 
over any cycles on the Schottky double. Now we are going to introduce new variables, 
IIq,, which can be called virtual a-periods of the differential S(z)dz on the Schottky dou- 
ble, since in all the cases when the Schwarz function has a meromorphic extension to 
the double they indeed coincide with the a-periods of the corresponding differential (see 
below in this section). 

Let us consider the differential 

cM k = dnt-J2 B k (z a )dW a (k > 1) ( 3.27) 

a 

where 

B k {z) = f dB k (3.28) 
jo 

is a polynomial of degree k. It is a meromorphic differential on £ with the only pole at 
oo on the upper sheet, where it has the form 

dn k (z) = dB k (z) + 0(z- 2 )dz (3.29) 

From (j3.21j) it is clear that the differential dQ k has vanishing a-periods 

/ dn k = (3.30) 

i.e. it is a canonically normalized meromorphic differential. The normal displacements 
of the boundary given by real and imaginary parts of the normal derivative d n Q k define 
a complex tangent vector field 

dl=dt k -Y.B k {z a )di ( 3. 31) 

a 

to the space of multiply-connected domains. These vector fields keep fixed the formal 
variable 

n a = (p a + 2ReJ2 B k (z a )r k (3.32) 
k 

In general situation this variable is only a formal one because the sum generally does not 
converge. Thus, we call n Q the virtual a-period of the Schwarz differential S(z)dz, since 
in the case when the Schwarz function has a meromorphic extension to the double S, the 
sum does converge and the corresponding quantity does coincide with the a a -period of 
the extension of the Schwarz differential. 



Elementary deformations and the operator V(z). Like in sect. 2, we introduce 
the elementary deformations 

5 a with 6n(0 = -id n G(a,£), a G D c 

1 (3.33) 
6^ with 8n{£) = - -d n u a {£) , a = l,...,g 
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where u a (z) is the harmonic measure of the boundary component b Q (see (I3.1J1 ). The 
deformations 6^ were considered in jTU] in connection with so called quadrature domains 

H2H2DI- 

In complete analogy with sect. 2 one can derive the following formulas for variations 
of the local coordinates under elementary deformations: 

2 , 5^r k = 0, 5^<f)p = -2e5 a(3 (3.34) 

The first two formulas are particular cases of 

6 a f h{QdX= e -l h(Od n G(a,()\d(\ = -enh(a) 
Jd" 2 JdD c 

which is valid for any harmonic function h in D c (the last equality is just the formula for 
solution of the Dirichlet problem). Similarly, 

6™ f h(()d 2 ( = e -l hd n u a \d(\ = e -i d n h \d(\ = -ie I d c hd( 

JD C 2 JdD' 2 Jb a JdD a 

(the Green formula was used), and the last two formulas in ()3.34j) correspond to the 
particular choices of h(z). Variations of the variables n a (in the case when they are 
well-defined) then read: 

5 z U a = 0, S^ILf, = -2e5 af3 (3.35) 

Therefore, for any functional X on the space of the multiply-connected domains the 
following equations hold: 

5 z X = eV{z)X (3.36) 
8 {a) x = -2edtX = -2ed^X (3.37) 
The differential operator V(z) in the multiply-connected case is defined by the formula 

v W = a£ + E {M*K + Mz)d? k ) (3 . 38) 
k>i 

The functional X can be regarded as a function X = X^((f) a , r^) on the space of the local 
coordinates 4> a , T k, or as a function X = X u (Jl a , r k ) on the space of the local coordinates 
n a , r fc . We would like to stress once again, that although in the latter case the variables 
U a are formal their variations under elementary deformations and the vector- fields , 
which keep them fixed, are well-defined. 

For completeness, let us characterize elementary deformations 5 a in terms of mero- 
morphic differentials on the Schottky double (as we have already seen, deformations 5^ 
correspond to holomorphic differentials). 

• Let d t { a ) be the vector field in the space of multiply-connected domains corresponding 
to the elementary deformation S a . Then the differential —(d t ( a )(d( — d t ( a )(d() ex- 
tends from dD c to the Schottky double S. Its extension is a meromorphic third-kind 
Abelian differential dVt^ a,a ^ which has simple poles at the points a and a on the two 
sheets of the Schottky double (with residues ±1) and vanishing ^-periods. 



5 a r k = eA k (a), 5 a <j> a = elog 



I — 
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In terms of the Green function we have: 

{2d z G(a, z)dz on the upper sheet 
—2dzG(a, z)dz on the lower sheet 
(cf. ()3.23|) and ()3.33j0 . Note that the differential dQ Q introduced before coincides with 

Let K(z, ()d( be a unique meromorphic Abelian differential of the third kind on S 
with simple poles at z and oo on the upper sheet with residues ±1 normalized to zero 
a-periods. (Note that as a function of the variable z it is multi- valued on £.) Then 

2d ( G(z, ()d( - 2«9 c G(oo, ()d( = K(z, ()d( + K(z, ()d( (3.39) 

and the differential dQ k {C) can be represented in the form 

dn k (() = ^<f K(u, QdB k {u) (3.40) 

where the w-integration goes along a big circle around infinity. Using (J3.14)) we obtain 
that 

- £ A k (z)dV k (() = £ I K{uX)du = K(z, QdC (3.41) 

Therefore, the following expansion of the derivative of the Green function holds: 

2d c G(z,()d( = dQ (C) - £ (A k (z)<m k (Q + A k Jz)dCl k (C)) (3.42) 

fc>i 

Here dfl k is a unique meromorphic differential on £ with the only pole at infinity on 
the lower sheet with the principal part —dB k (z) and vanishing a-periods. This formula 
generalizes eq. (3.8) from ,2] to the multiply-connected case. 



The F-function. Applying the variational formulas (jl.4j) . (j3.3|) . ()3.4|) . we can find 
variations of the Green function, harmonic measure and period matrix under the elemen- 
tary deformations. In this way we obtain a number of important relations which connect 
elementary deformations of these objects: 

S a G(b, c) = SbG(c, a) = 5 c G(a, b) 
5 a u a (b) = 5^G(a,b) = 5 b u a (a) 

6Ww p (z) = 5^u a (z) (3.43) 
5 z Q aP = 8^u p {z) 

8^n M = s^n ja = 5^n a(3 
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From ()3.36|) . ()3.37|) it follows that the formulas ()3.43|) can be rewritten in terms of the 
differential operators V(z) and d a := d/d(f> a = d/dH a : 

V{a)G{b,c) = V(b)G{c,a) = V{c)G{a,b) 
V{a)u a (b) = -2d a G(a,b) 

d a up(z) = d p u) a (z) (3.44) 
V{z)Vt aP = -2d a up(z) 

These integrability relations generalize formulas (j2.15|) to the multiply-connected case. 
The first line just coincides with f!2.15|) while the other ones extend the symmetry of the 
derivatives to the harmonic measure and the period matrix. 

Again, ()3.44|) can be regarded as a set of compatibility conditions of an infinite hi- 
erarchy of differential equations. They imply that there exists a function F = F(U a , r) 
such that 

G(a, b) = log (a" 1 - 6 _1 | + ^V(a)V(6)F (3.45) 
Ua ( z ) = -d a V(z)F (3.46) 
T aP = M aP = 2m d a d p F (3.47) 
The function F is the (logarithm of the) tau-function of multiply-connected domains. 



Dual moments and integral representation of the tau-function. To obtain the 
integral representation of the function F, we proceed exactly in the same manner as in 
sect. El (see also [2 j for more details). 

Again, set <t(z) = V(z)F. Eqs. (13.45)) and (|3.46|) determine the function &(z) for 
z G D c via its variations under the elementary deformations: 



6 a $(z) = -2e\og 



2eG(a, z) 



(3.48) 



It is easy to verify that the function 

2 



$(z) = / \og\z- 1 - C l \d 2 C = ~ $(0) + r log |^| 

7T JD 

satisfies ([3.48)1 . Indeed, using ()3.33j) . variation of ()3.49j) reads 



/ io g |^- i -r 1 |rf 2 c) = 1 / i^ia^oiogi^-r 1 ! 

JD / 7T JdD" 



e 

7T J dD 



(3.49) 



/ \dnd n G(a,0 (\og\z- 1 -C 1 ]- G(z,0) = -2e\og a' 1 -z- 1 +2eG(a,z 
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where we have used properties of the Dirichlet Green function and the fact that Dirichlet 
formula restores harmonic function from its value at the boundary. Similarly, for z G D c 
we obtain: 



■ i -- / \og\z- i -c x \sc] = -i md n u a (ti)\og\z~ i - 1 

JD J 71 JdD c 



-1| 



71 JD 

e 

71 JdD 



I \d£\d n u a (£) (logl^-T 1 ! - G(z,0 

JdD c v 

= - / \d£\d n (log {z- 1 - T 1 ! - G(z, 0) = 2eu a (z) 

71 Jb a V 7 

The same calculation for z G D yields 

f if z G D 
= J (3.50) 

[ 2e5 a(3 if z G D^, /3 = 

We see that expression in ()3.49J) coincides with $ given by (|2.32() . where D is now 
understood as the union of all D a 's. The coefficients of an expansion of <3> at infinity 
define the dual moments v k : 

V(z)F = $(z) = -- [ \og\z- 1 -C 1 \d 2 ( = vo + 2ReY,VkA k (z) ( 3 .51) 

The coefficients in the r.h.s. of (|3.5H1 are moments of the union of the interior domains 
with respect to the dual basis 

v k = -l B k {z)d 2 z (3.52) 



7T JD 

From equation (|3.51|) it follows that 

v k = S%F (3.53) 

The same arguments show that the derivatives 

s a :=-d a F (3.54) 

are just areas of the holes (|3.11|) . Indeed, eqs. (|3.46|) . (|3.47|) determine these quantities 

via their variations: 5 a s a = euj a (a), S^s a = eVt a p. A direct check, using ()3.33|) . shows 



that 

'1 1 



d , c daF = _ MM (3.55) 

7T 



V(z)F = — [ log 

7T JD Z Q 

For example, 

smf = ^ ( [ $(z)d 2 z) = ~l md n u a (om + ^ / * w $(o d\ 

ATX \JD J 471 JdD" 171 JD 

In the last term we use (13.501) and obtain the result: 



S&F =-^-i \d£\dj>(£) + - f d 2 ( = ~ ( A$ d 2 C + es a = 2es a 

ATT Jb a 7T JD a ATT JD a 

(here A = 4<9 Z <9 2 is the Laplace operator). 

The integral representation of F is found in the same way through its variations which 
are read from (|3.55|) . The result is given by the same formulas (12.35(1 and (|2.36(1 as in 
the simply-connected case with the understanding that D = U^ =0 D a is now the union of 
all D 's. 
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Algebraic domains. In what follows we restrict our analysis by the class of algebraic 
domains. In the simply-connected case dealt with in the previous section the algebraic 
domains are simply images of the exterior of the unit disk under one-to-one conformal 
maps given by rational functions whose singularities are all in the other "half" of the 
plane, i.e. inside the unit circle. Note that the boundary of the unit circle is the set 
of fixed points of the inversion w — > 1/w which is the antiholomorphic involution of the 
w-plane compactified by a point at infinity (the Riemann sphere) . 

Planar multiply-connected algebraic domains can be defined as the domains for which 
the Schwartz function has a meromorphic extension to a higher genus Riemann surface 
(a complex algebraic curve) with antiholomorphic involution. More precisely, let E be 
a real Riemann surface by which we mean a complex algebraic curve of genus g with 
an antiholomorphic involution such that the set of fixed points consists of exactly g + 1 
closed contours (such curves are sometimes called M-curves). Then £ can be naturally 
divided in two "halves" (say upper and lower sheets) which are interchanged by the 
involution. Algebraic domains with g holes in the plane can be defined as images of 
the upper half of the real Riemann surface under bijective conformal maps given by 
rational (meromorphic) functions on E (see, e.g. |21| ). For the purpose of this paper, it 
is convenient to use another, more direct characterization of algebraic domains. 

The domain D c is algebraic if and only if the Cauchy integrals (j3.15j) 

C a {z) = ±-J forzeD a 

ZTTt JdD Q — Z 

are extendable to a rational (meromorphic) function J(z) on the whole complex plane 
with a marked point at infinity (see [21] )• R is important to stress that this function is 
required to be the same for all a. The equality 

S( z ) = J(z)-C-(z) 

valid by definition for z G <9D C can be used for analytic extension of the Schwarz function. 
The function C~{z) is analytic in D c . Therefore, J(z) and S(z) have the same singular 
parts at their poles in D c . One may treat S(z) as a function on the Schottky double 
extending it to the lower sheet as z. 

It is also convenient to introduce 

Viz) = f J(z)dz (3.56) 
Jo 

which is multi-valued if J(z) has simple poles (to fix a single-valued branch, we make 
cuts from oo to all simple poles of J(z)). In fact we need only the real part of V(z). In 
neighborhoods of the points z a one has 

J(z)dz = J2 i~kdB k (z), V(z) = T kB k (z) /g 57 n 

k>l k>l 

The formula (EP2D 

U a = a + 2ReV(z a ), z a G D a (3.58) 

shows that for the algebraic domains the variables n a , introduced in the general case as 
formal quantities, are well-defined. It is easy to show that they are equal to the a-periods 



26 



of the differential S(z)dz on the Schottky double E. Indeed, using the fact that C°(z) 
and C a (z) represent restrictions of the same function J(z), one can rewrite (|3.18j) in the 
form 

(f) a = -2 Re (j* J(z)dz + J^{z + C~(z) - J(z))dz^j 

Under the second integral we recognize the Schwarz function. Combining this equality 
with the definition of U a ()3.32|) . we obtain: 



n Q = 2 Re [°(S(z) - z)dz = [° (S(z)dz - zdS(z)) = I S{z)dz (3.59) 



As an example of algebraic domains, it is instructive to consider the case when only a 
finite-number of the moments r k are non-zero. Let An be the space of multiply-connected 
domains such that 

r k = 0, k > N (3.60) 
Then the arguments similar to the ones used above show that 

• S(z) extends to a meromorphic function on E with a pole of order N — 1 at oo and 
a simple pole at ob. 



The function z extended to the lower sheet of the Schottky double as S(z) has a simple 
pole at oo and a pole of order N — 1 at ob. For a domain D c G An the moments with 
respect to the Laurent basis (cf. (j2.1|) ) 

t k = -—r I z- k d 2 z (3.61) 



coincide with the coefficients of the expansion of the Schwarz function near oo: 

N 



S(z) = Y J kt k z k - 1 + 0(^ 1 ), z^oo (3.62) 



k=l 



The normal displacement of the bondary of an algebraic domain, which changes the 
variable t k keeping all the other moments (and n Q ) fixed is defined by normal derivative 
of the function 2 Re f z dfl k . Here dQ k is a unique normalized meromorphic differential 
on E with the only pole at oo of the form 

dQ k = d(z k + 0(z' 1 )), I dQ k = (3.63) 

J 3 a 

Note that the differential dQ k is well-defined for a generic, not necessarily algebraic 
domain. Therefore, the normal derivative of the function 2 Re f z dfl k defines a tangent 
vector field on the whole space of multiply-connected domains. 

The space An is a particular case of algebraic orbits of the universal Whitham hier- 
archy. In this case the general formula (7.42) from [3] for the r-function of the Whitham 
hierarchy, after proper change of the notation, acquires the form 

1 1 9 

2F = --r 2 + t v + - ]T(2 - k)(r k u k + r k v k ) - £ U a s a (3.64) 

1 1 k>l a=l 
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which is a quasi-homogeneity condition obeyed by F (compare with formula (5.11) from 

El)- 

Let dQo be a unique normalized meromorphic differential on £ with simple poles at 
the infinities oo and 6b. Its Abelian integral 

\ogw(z) = [ dVl (3.65) 

defines in the neighborhood of oo a function w(z) which has a simple pole at infinity. The 
dependence of the inverse function z(w) on the variables t k is described by the Whitham 
equations for the two-dimensional Toda lattice hierarchy. These equations have the form 

d*z(w) = {n k (w), z(w)} := ^^d to zH - dMw)^ (3-66) 

Algebraic domains of a more general form correspond to the universal Whitham hierarchy. 
Let An 1 ,...,n 1 be the space of domains such that the extension of the Schwarz function 
S(z) to D c has poles of orders Nj — 1 at some points Zj (which possibly include oo and 
6b). Then, according to 0j, the variables 

r z j 1 

to,j = / S(z)dz, tkj = —res Zj (z — Zj) k ^ 1 S(z)dz 1 k — 1, . . . , Nj — 1 (3.67) 
J Co k 

together with the variables s a (or IT a ) provide a set of local Whitham coordinates on the 
space An 1 ,...,Ni- Note that the definition of the algebraic orbits of the universal Whitham 
hierarchy is a bit more general than the definition of algebraic domains given above. It 
corresponds to the case when the differential dS of the Schwarz function is extendable 
to D c as a meromorphic differential (in such domains are called Abelian domains). 
For example, let A^ be the space of multiply-connected domains such that 

T fe {1) = <f A k dS = I S\z)A k dz = -<f A' k S{z)dz = 0, k > N (3.68) 



/<9D JdD J dD 

This spaces is characterized by the following property: there are constants K a such that 
S(z) + K a extends to a meromorphic function on D c with a pole of order iV at oo. The 
variables ^ 

K a , s a , tk = -—r f z~ k zdz, k > 1, (3.69) 

llilk JdD 

are local coordinates on A^ . 

The two cases when S(z) or its derivative S'(z) have a meromorphic extension to 
D c are particular examples of the whole hierarchy of integrable domains, which can be 
defined in a similar way by the condition that the m-th order derivative of the Schwarz 
function S(z) admits a meromorphic extension to the Schottky double E. For example, 

(2) 

let An be the space of multiply-connected domains such that 

T, (2) = / S"(z)A k dz = -l A'lS(z)dz = 0, k > N (3.70) 

J 3D JdD 

This space is characterized by the following property: there are linear functions k a (z) = 
K® + K\z such that S(z) + k a (z) extends to a meromorphic function on D c with a pole 
of order N + 1 at oo. The variables 

K a (z), s a , t k = — — - I z~ k zdz, k>l, (3.71) 
2mk JdD 
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are local coordinates on A N . The other spaces A^ 1 with m > 2 can be defined in a 
similar way. 



4 The duality transformation 

The independent variables II a ()3.59|) or <p a used in the previous section are not as trans- 
parent as the dual variables s a (j3.11j) . which are simply areas of the holes D^. In this 
section we show how to pass to the set of independent variables si,...,s g (|3.11|) (together 
with the infinite set of Tfc's). This transformation is similar to the passing from "external" 
to "internal" moments in the simply-connected case (see sect. 5 of 2J). The difference 
is that only a finite number of times are subject to the transformation while the infinite 
set of Tfc's remains the same. 

The change to the variables s a can be done in general case of domains with smooth 
boundaries. However, it is the change Il a — > s a rather than 4> a — > s a that leads to a 
transparent duality. Since n a 's are only defined as formal ("virtual") variables for do- 
mains with smooth boundaries, we shall restrict our consideration to the class of algebraic 
domains discussed at the end of the previous section. In this case the variables U a are 
well defined. 



The Legendre transform. Passing from IT a to s a is a particular duality transforma- 
tion which is equivalent to the interchanging of the a and b cycles on the Schottky double 
S. This is achieved by the (partial) Legendre transform F(H a , r) — > F(s a , r), where 

F = F+]Tn Q s Q (4.1) 

a=l 

The function F is the "dual" tau-function. Below in this section, it is shown that F 
solves the modified Dirichlet problem and can be identified with the free energy of a 
matrix model in the planar large N limit in the case when the support of eigenvalues 
consists of a few disconnected domains (a so called multi-support solution, see [22] and 
references therein). 

The main properties of F follow from those of F. According to (|2.34j) . (J3.55j) we have 
dF = — Y^ a s a dH a + Ysk Vkdrk (for brevity, k is assumed to run over all integer values, 
T_fc = ffe, etc.), so dF = J2 a n a ds a + J2k Vkdru. This gives the first order derivatives: 

9F dF , A , 

n Q = j— , v k = — 4.2 
ds a dr k 

The second order derivatives are transformed as follows (see e.g. [23J). Set 

d 2 F „ d 2 F „ 8 2 F 



al3 dliadlip ' ah dU a dr k ' lk dridr k 



29 



and similarly for F. Then 

Fa/3 = — (F -1 )^ 



v ak ^ ' (F ^)ot~fF^k 
7=1 



(4.3) 



Fik Fife ^ ' Fi^i^F )t7' Fly'k 
7,7'=1 

Here (-F _1 ) Q /3 means the matrix element of the matrix inverse to the g x g matrix F a p. 

Using these formulas, it is easy to see that the main properties (13.45)) . ()3.46|) and 
(13.47)) of the tau-function are translated to the dual tau-function as follows: 

G(a, b) = log (a" 1 - ft -1 1 + ^V(a)V(o)F (4.4) 

2mu a {z) = -d s y(z)F (4.5) 
d 2 F 

2mT aP =-—— (4.6) 

0S a 0Sf3 

where r^-derivatives in V(z) are taken at fixed s a . The objects in the left hand sides of 
these relations are: 

G(a,b) = G(a,b) + in Y u a (a)f a pu)p(b) (4.7) 

9 

^a(z) = T aP u) P {z) (4.8) 
f = -T- 1 = -IT 1 (4.9) 

7T 

The function G is the Green function of the modified Dirichlet problem to be discussed 
below. The matrix T is the matrix of a-periods of the holomorphic differentials dW a on 
the double S (so that u) a (z) = W a (z) + W a (z)), normalized with respect to the b-cycles 

-I dW p = 8 a p, I dWp = f aP 

i.e. more precisely, the change of cycles is a a — > b a , b Q — > -a a . 

An important remark is in order. By a simple rescaling of the independent variables 
one is able to write the group of relations (J4.4j) - (J4.6)) for the function F in exactly the 
same form as the ones for the function F (13.45)1 — (|3.47jl . so that they differ merely by the 
notation. We use this fact in sect. El 
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The modified Dirichlet problem. The modified Green function ()4.7|) solves the 
modified Dirichlet problem which can be formulated in the following way. 

One may eliminate all except for one of the periods of the Green function G, thus 
making it similar, in this respect, to the Green function of a simply-connected domain 
(recall that the latter has the non-zero period 27r over the only boundary curve b ). This 
leads to the following modified Dirichlet problem (see e.g. [21] ) : given a function Uo(z) 
on the boundary, to find a harmonic function u(z) in D c such that it is continuous up 
to the boundary and equals uo(z) + C a on the a's boundary component. Here, C a 's are 
some constants. It is important to stress that they are not given a priori but have to 
be determined from the condition that the solution u(z) has vanishing periods over the 
boundaries b>i, . . . , b g . One of these constants can be put equal to zero. We set Cq = 0. 
This problem also has a unique solution. It is given by the same formula (jl.lj) in terms 
of the modified Green function G(z,Q. The definition of the latter is similar to that of 
the G(z, () but differs in two respects: 

(Gl) G(z, () is required to have zero periods over the boundaries bi, . . . , b s ; 

(G2) The derivative of G(z, () along the boundary (not G(z, () itself!) vanishes on the 
boundary. 

Under the condition that G(z, () = on bo such a function is unique. The function given 
by (j4.7j) just meets these requirements. We conclude that the modified Green function is 
expressed through the dual tau-function F. 

Note that variations of the modified Green function under small deformations of the 
domain are described by the same Hadamard formula (jl.4j) . where each Green function is 
replaced by G. This follows, after some algebra, from the formula for G in terms of G, u a 
and Q a /3- Therefore, all the arguments of sect.EJcould be repeated in a completely parallel 
way starting from the modified Dirichlet problem. One may also say that the functions G 
and G differ merely by a prefered basis of cycles on the double: the differential d z Gdz has 
vanishing periods over the a-cycles while d z Gdz has vanishing periods over the b-cycles. 

The function w a (z) such that G(z, a) = — log |u> a (;2)| maps D c onto the exterior of the 
unit circle which is slit along g concentric circular arcs (see fig. El)- Since the periods of 
the function G vanish, the function w a is single-valued. Positions of the arcs depend on 
the shape of D c as well as on the point a which is mapped to oo. The radii of the arcs, 
R a , are expressed through the dual tau-function as 

log R 2 a = d s y(a)F (4.10) 

In particular, for a = oo we have logi?^, = d Sa d TQ F (cf. eq. ()2.20|) for the conformal 
radius). 



Relation to multi-support solutions of matrix models. The partition function 
of the model 3 , written as an integral over eigenvalues, reads: 



<N 



1 r ( N 1 N \ N 

- / exp £log {Zi - + - £ U{ Zi ) I] Sz 3 (4.11) 



One may have in mind the model of all complex or mutually hermitean conjugated or normal 
matrices. 



31 




Figure 5: The image of a triply-connected domain under the conformal map w a (z) such 
that G(z, a) = — log |w a (z)|. The function w a maps the domain onto the exterior of the 
unit (dashed) circle with g = 3 concentric circular cuts. Positions and lengths of the arcs 
depend on the shape of D c in fig. 01 and depend also on the normalization point z = a 
which is mapped to oo. 



The matrix model potential U is usually chosen to be of the form 



U{z) = -zz + V{z) + V(z) (4.12) 

where V(z) is at the moment some polynomial; however, we will see immediately that 
the coincidence of the notation with ()3.57|) is not accidental. The parameter H, in the 
large N limit, tends to zero simultaneously with iV — > oo in such a way that to = Nh is 
kept finite and fixed. In the leading order, one can apply the saddle point method. The 
saddle point condition is 

N 1 



The density of eigenvalues is a sum of two-dimensional delta-functions: 

p(z)=nhY / ^ 2 \z-z l ) 

i 

In the large N limit, one treats it as a continuous function normalized as - / p(z)d 2 z = t . 
In terms of the continuous density the saddle point equation reads 

1 f p{Od2(: + d z U(z) = (4.13) 



7C J Z-( 

The solution is well known and easy to obtain: the extremal density po(z) is constant in 
some domain D (the support of eigenvalues which can be a union of disconnected domains 
Dq,) and zero otherwise. More precisely, 



Po{z) 



l z e D 

z6C\DeD c 



Note that the saddle point condition is imposed only for z inside the support of eigen- 
values. Writing (cf. (13371) ) 

V{z) = j2nB k {z) 
k=i 
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multiplying eq. ()4.13|) by A k {z) and integrating it over the boundary of the support of 
eigenvalues, one finds that the domain D is such that the coefficients r k are moments 
of its complement with respect to the basis functions A k and the higher moments (with 
numbers greater than p) vanish. As is proven in sect. El these conditions, together with 
the normalization condition for the density, locally determine the shape of the support of 
eigenvalues. Equivalently, one might parametrize the polynomial as V(z) = J2k=i^kZ k , 
then tk are moments of the C \ D with respect to the functions z~ k . 

An interesting problem is to obtain, for a given value of to, necessary and sufficient 
conditions on the polynomial V for the support of eigenvalues to be a union of g + 1 
disconnected domains ("droplets") with non-zero filling. One may approach this problem 
from a "classical" limit of very small (point-like) droplets. Clearly, for our choice of the 
signs in (|4.11|) and (j4.12j) . the stable point-like droplets are located at minima of —U(z), 
or equivalently at maxima of U(z). As soon as we use the basis which explicitly depends 
on the marked points z a G D a , it is natural to consider the germ configuration with point- 
like eigenvalue droplets at the points z a . It is easy to see that the sufficient conditions 
for the potential (|4.12j) to have maxima at the points 

z a = V'(z a ), \V"(z a )\ < 1 for all a (4.14) 

The first one means that there is an extremum of the potential U(z) at the point z a . The 
second one ensures that eigenvalues of the matrix of second derivatives of the potential at 
the extremum are both negative, i.e. the extremum is actually a maximum of U(z). The 
first condition literally coincides with the one used in [22] for the completely degenerate 
curve (point-like droplets) while the second one requires that not all extrema are filled, so 
that the "smooth" genus g is always less than the maximal possible genus (p— l) 2 — 1. Let 
now Vq(z) be a minimal degree polynomial that obeys these conditions. Then perturbed 
polynomials of the form 

V(z) = V (z)+ ]T r k B k (z) 

k>g+2 

obey the same conditions for sufficiently small r*., and this is the advantage of using the 
basis (J3.1(J|) . (J3.13|) from the perspective of matrix models. 

The saddle point equation (|4.13|) just means that the "effective potential" - J D log \ z — 
C\ 2 d 2 ( + U(z) is constant in each D a , i.e. for z G D a it holds: 

- / log \z - CM 2 C + U(z) =v + U a (4.15) 

7T JD 

where Vq = - J* D log \z\d 2 z, and so in our normalization Ho = 0. 

Let N a be the number of eigenvalues in D a , then limTv^oo^^a — s a . First we find 
the free energy with some fixed s a : 

= hl I Po{z) lo gl^-C|Po(C)^ 2 C + ^/ Po (z)U(z)d 2 z 
where one should substitute (|4.15jl . The result is 

= -—J f log\z- 1 -C 1 \d 2 zd 2 ( + Y / Il a s a = F(s a ,r) (4.16) 

s a fixed 7T JDJD 



lim (h 2 \ogZ N 



lim \h 2 \ogZN 

N— >oo \ 
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This quantity depends on s a 's. It is given by the value of the integrand in (j4.11j) at the 
saddle point with fixed s a . 

If one wants to take into account the "tunneling" of eigenvalues between different 
components of the support, s a are no longer free parameters, and the free energy of 
the "planar limit" of the matrix model should be obtained by extremizing F with respect 
to s a (with fixed t ). It follows from the above that d Sa F = U a and so the extremum is 
at U a = 0. Therefore, 



(4.17) 



F<® = 

n Q =o 

where F is given by (|2.35|) . The results of sect. El imply that the growth of the support of 
eigenvalues, when N —> N + 5N and the tunneling is taken into account, is given by the 
normal displacement of the boundary proportional to the normal derivative of the Green 
function with the minus sign, like in (|2.9J1 . Since this quantity is always non-negative, we 
conclude that if a point belongs to the support of eigenvalues, it does so as t increases. 
In particular, if one starts with point-like droplets at the points discussed above, 

then these points always remain inside the droplets. 

Let us note that different aspects of multi-support solutions of the 2-matrix model 
and matrix models with complex eigenvalues were discussed in [2"2" | I26 | l2"H]. 

For matrix models one usually restricts oneself to the algebraic case of finite amount 
of nonvanishing moments playing the role of the coefficients of the matrix model po- 
tential (|4.12|) . In such case the corresponding complex curve can be described by an 
algebraic equation 22 j, which can be thought of as auxiliary constraint to the second 
derivatives of F. These auxiliary constraints look similar to the reduction conditions in 
the case of Landau- Ginzburg topological theories (see, for example, discussion of such 
conditions in the context of dispersionless Hirota equations and WDVV equations in [7j). 
Their meaning is that the derivatives w.r.t. the times {r^} for k > g + 1 (where g is 
genus of corresponding algebraic complex curve) can be expressed through the derivatives 
restricted to k < g + 1. 



5 Green function on the Schottky double and gen- 
eralized Hirota equations 

Let us now turn to the generalization of the dispersionless Hirota equations to the 
multiply-connected case. 

For a unified treatment of the two "dual" representations of the Dirichlet problem 
discussed in the two previous sections, we make a simple change of variables. Namely, 
let us introduce the generic "period" variables X a which are identified either with n a or 
2m s a depending on the choice of the set of cycles, and the function J-(X a , t) equal to 
F(X a ,r) or F(^,r) respectively. Then the main relations (|B~43J) - (ETT7I) and (jOjl - 
(14. 6|) acquire the form 

G(a, b) = log (a" 1 - 6 _1 | + ~V(a)V(fc);F (5.1) 
uj a (z) = -9 Q V(^)J r (5.2) 
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T af) = 2irid a dpJ 7 (5.3) 

where d a := dx a and G, u> and T stand for the corresponding objects with or without 
tilde, depending on the chosen basis of cycles. 

We will see below that, in analogy to the simply-conected case, any second order 
derivative of the function T w.r.t. r k (and f k ), jF ifc , will be expressed through the deriva- 
tives {J-'ap} where a, (3 = 0, . . . ,g together with {.F QTi } and their complex conjugated. 
To be more precise, one can consider all second derivatives as functions of {J 7 ^, J 7 ^} 
modulo certain relations on the latter, like the relation (j5.14|) to be discussed below. 
Sometimes on this "small phase space" more extra constraints arise, which can be writ- 
ten in the form similar to the Hirota or WDVV equations |2Z] ; we are not going to discuss 
this issue here, restricting ourselves to the generic situation. 



The Abel map. To derive equations for the function F = T in sect. EJ we used the 
representation ()2.23j) for the conformal map w(z) in terms of T and eq. (jl.2)l relating 
the conformal map to the Green function which, in its turn, is expressed through the 
second order derivatives of T. In the multiply-connected case, our strategy is basically 
the same, with the suitable analog of the conformal map w(z) (or rather of log w(z)) 
being the embedding of D c into the g- dimensional complex torus Jac, the Jacobi variety 
of the Schottky double. This embedding is given, up to an overall shift in Jac, by the 
Abel map z i— > W(z) := (Wi(z), . . . , W g (z)) where 




is the holomorphic part of the harmonic measure u a . By virtue of (|5.2|) . the Abel map 
is represented through the second order derivatives of the function T: 

W a (z) - W a (oo) = [ Z dW a = -d a D(z)F (5.5) 

J oo 

2ReW„(oo) = co> Q (oo) = — <9 TO <9 Q .F (5.6) 
The last formula immediately follows from (|5.2j) . 



The Green function and the prime form. The Green function of the Dirichlet 
boundary problem, appearing in (|5.1)) . can be written in terms of the prime form (jA.4|) 
on the Schottky double (cf. JOJ)): 



G(z,C)=log 



E(z,C) 



E(z,C) 



(5.7) 



Here by ( we mean the (holomorphic) coordinate of the "mirror" point on the Schottky 
double, i.e. the "mirror" of £ under the anti holomo rphic involution. The pairs of such 
mirror points satisfy the condition W a (() + W a {Q = in the Jacobian (i.e., the sum 
should be zero modulo the lattice of periods). The prime form 4 is written through the 

4 Given a Riemann surface with local coordinates 1/z and 1/z we trivialize the bundle of —\- 
diffcrcntials and "redefine" the prime form E(z, Q — + E(z, ()(dz) 1 / 2 (dQ 1 ^ 2 so that it becomes a function. 
However for different coordinate patches (the "upper" and "lower" sheets of the Schottky double) one 
gets different functions, see, for example, formulas (|5.8|) and (|5.9|) below. 
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Riemann theta functions and the Abel map as follows: 



E{ZX) ~ h(z)h(0 (5 ' 



when the both points are on the upper sheet and 



, Ts 6JW(z)+W(()) , s 

= — — - 5.9 

i/»(z)MO 

when 2! is on the upper sheet and ( is on the lower one (for other cases we define E(z, Q = 
E(z, C), E(z,C) = E(z,Q). Here 0,(W) = 6g*(W\T) is the Riemann theta function 
flA.2|) with the period matrix T a p = iTtidadpT and any odd characteristics 8*, and 

h\z) = -z 2 £ 8*A Q ) d *W a (z) = z 2 j2 0*, a (O) £ ^(z)^^ (5.10) 

o=l a=l k>l 

Note that in the l.h.s. of (J5.9|) the bar means the reflection in the double while in the 
r.h.s. the bar means complex conjugation. The notation is consistent since the local 
coordinate in the lower sheet is just the complex conjugate one. However, one should 
remember that E(z, Q is not obtained from ()5.8|) by a simple substitution of the complex 
conjugated argument. On different sheets so defined prime "form" E is represented by 
different functions. In our normalization (J5.9j) iE(z, z) is real (see also Appendix B) and 

lim EM -< 



In particular, lim^^oo zE(z, oo) = 1. 



The prime form and the tau-function. In (|5.7|h the /i-functions in the prime forms 
cancel, so the analog of (|2.17jl reads 



log 



0,(W(z)-W(O) 



^(W(z)+W(C)) 



2 


1 1 


= log 









V(z)V(0^ 



(5.11) 



This equation already explains the claim made in the beginning of this section. Indeed, 
the r.h.s. is the generating function for the derivatives Tik while the l.h.s. is expressed 
through derivatives of the form T a h and T a p only. The expansion in powers of z, £ allows 
one to express the former through the latter. 

The analogs of eqs. (|2.19|) . ()2.20|) are, respectively: 

2 



log 



0«,(W(z) -W(oo)) 



ejw(z) 



log 



W(oo)) 
h 2 (oo 



-\og\z\ 2 + d T0 V(z)J r 



Here uj(z) 



2ReW(^) = 
h 2 (oo) 



0.(«(oo)) 
(uJi(z), . . .,U) g (z)) and 



d 2 T 

U T J 



(5.12) 
(5.13) 



lim z6* 

z— >oo 



[ K dw) = - y £e.< a (o)d a d n r 

Jo ° ' a=l 
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A simple check shows that the l.h.s. of ()5.13|) can be written as —2log(iE(oo,6b)). As 
is seen from the expansion G(z, oo) = — log|z| — log(iE(oo, do)) + O^z^ 1 ) as z — ► oo, 
(iE(oo, do)) -1 is a natural analog of the conformal radius, and (|5.13j) indeed turns to 
()2.20|) in the simply-connected case (see Appendix B for an explicit illustrative example) . 
However, now it provides a nontrivial relation on T a ps and T a iS: 



Y.B^dad^ (y^e^dpd^A = ^(o;(oo))e^ 



(5.14) 



so that the "small phase space" contains the derivatives modulo this relation. 

The next steps are exactly the same as in sect. |2J we are going to decompose these 
equalities into holomorphic and antiholomorphic parts. The results are conveniently 
written in terms of the prime form. The counterpart of (|2.22jl is 

Tending ( — > oo, we get: 

log ^. Z,00 \ =logz + log£(oo, oo) - d T0 D(z)F (5.16) 
E(z, oo) 

Separating holomorphic and antiholomorphic parts of ()5.15j) in (, we get analogs of ()2.24j) 
and 42251): 

-■og |fefe^ = gM 5( C > (5.18) 
E{z, oo)£/(oo, C) 

Combining these equalities (with merging points z —> ( in particular), one is able to 
obtain the following representations of the prime form itself: 

E(z, C) = (z- 1 - C -i) e -K^)-s(C))^ (5.19) 

%E{z, C) = e -K^+^(^)+5(C))^ (5.20) 

Note also the nice formula 

iE(z,z) = e -5 v2 (^ (5.21) 



Generalized Hirota relations. For higher genus Riemann surfaces there are no simple 
universal relations connecting values of prime forms at different points, which, via (J5.19)) . 
(I5.2(J|) . could be used to generate equations on T. The best available relation is the 
celebrated Fay identity (jA.5|) . Although it contains not only prime forms but Riemann 
theta functions themselves, it is really a source of closed equations on J 7 , since all the 
ingredients are in fact representable in terms of second order derivatives of T in different 
variables. 

An analog of the KP version of the Hirota equation ()2.26|) for the function T can 
be obtained by plugging eqs. (|5.5j) and (|5.19|) into the Fay identity (jA.5J) . As a result, 
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one obtains a closed equation which contains second order derivatives of the T only 
(recall that the period matrix in the theta-functions is essentially the matrix of the 
derivatives T a p). A few equivalent forms of this equation are available. First, shifting 
Z — > Z — W 3 + W 4 in (jA.5|) and putting z± = oo, one gets the relation 

zU 



(a - b)e D{a)D{b):F e[ [ a dw+ [ dW-Z) e ( f dW 

\Joo J oo J \J oo 

+ (b- c ) e D i b )D{c)r e ^J b dW + J c dW _ Z ^J o (^J a rfW-zj + (5.22) 
+ (c - a)e D{c)D{a)T 6 ( f dW + r dW -Z) 9 I f dW - Z J =0 

\J oo J oo / \Joo J 

The vector Z is arbitrary (in particular, zero). We see that ()2.26|) gets "dressed" by the 
theta-factors. Each theta-factor is expressed through T only. For example, 

6 ( [ Z dW) = ex P ( " 27r2 E n <*npdl p F - 2m £ n a d a D{z)T 

KJo ° 7 n Q eZ V af3 a j 

Another form of this equation, obtained from ()5.22|) for a particular choice of Z, reads 
(a-6) C - 1 e^ (D(a)+D(fe))2 ^(c)^ ( r dW+ C dW] + [cyclic per-s of a, 6, c] = (5.23) 

\J oo J oo I 

Taking the limit c — * oo in ()5.22|) . one gets an analog of ()2.27|) : 



D{a)D{b)F 

" e(/-rfw-z)e(/idw-z) 



(5.24) 



a-b T1 a-b^dZ a ed^W-Z) 

which also follows from another Fay identity ()A.6jl . 

Equations on T with ffc-derivatives follow from the general Fay identity (jA.5|) with 
some points on the lower sheet. Besides, many other equations can be derived as various 
combinations and specializations of the ones mentioned above. Altogether, they form 
an infinite hierarchy of consistent differential equations of a very complicated structure 
which deserves further investigation. The functions T corresponding to different choices 
of independent variables (i.e., to different bases in homology cycles on the Schottky 
double) provide different solutions to this hierarchy. 



Higher genus analogs of the dispersionless Toda equation. Let us show how the 
simplest equation of the hierarchy, the dispersionless Toda equation ()2.30|) . is modified 
in the multiply-connected case. Applying d z d^ to both sides of (j5.18|) and setting ( = z, 
we get: 

(dD(z))(8D(z))F = -d z d- z hg E(z,z) 

Here dD(z) is the z-derivative of the operator D(z): dD(z) = J2k A' k (z)d Tk . To transform 
the r.h.s., we use the identity (jA.8|) (Appendix A) and specialize it to the particular local 
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parameters on the two sheets: 



\z\ d z d s \og E(z,z) 



6 2 (Z)E 2 (z } z) 



+ \z\ A J2( l °S^)),apd z W a (z)d,W p (z) 

<x,P 



Tending z to oo, we obtain a family of equations (parametrized by an arbitrary vector 
Z) which generalize the dispesrionless Toda equation for the tau-function: 



d Tx d fx T 



e 2 (z) 



< ■ -J2( l °sO(Z)) tal3 (d a d T1 T)(d (} d fl T) (5.25) 

0,(9=1 



(Here we used the z — > oo limits of ()5.5|) and (|5.21jl .) The following two equations 
correspond to special choices of the vector Z: 



d Tl d fl F+ £ {\og6{V)) ia p{d a d T1 F){dpd fl T) = g2( ^ (00)) eV 



# 2 (0) 



(5.26) 



a,P=l 



(5.27) 



Finally, let us specify the equation (|5.25|) for the genus g — 1 case. In this case there is 
only one extra variable X := X\ (111 or 2ms\). The Riemann theta-function 9(u)(oo) + Z) 
is then replaced by the Jacobi theta-function i? (d x d To !F — Z\T) = $3 (dxd TQ J- — Z\ T), 
where the elliptic modular parameter is T = 2m d\J z ) and the vector Z = Z has only 
one component. The equation has the form: 



tf 3 (d x d T0 F+Z\ 2m d\T) tf 3 {d x d TQ T-Z\ 2m d\T) d , T 



d\(Z\2md\T) 
- d\\og^[z\2md\f) (d x d Tl F)(d x d fl F) 
Note also that equation (J5.14j) acquires the form 



e to* 



(5.28) 



{d x d n F)(d x d fl F) 



^(d x d T0 ^\2md 2 x ^y 
^' 1 (0\2md x J 7 ) 



a 2 .f 
e T o 



(5.29) 



where 1?* = $1 is the only odd Jacobi theta-function. Combining ()5.28|) and fj5.29j) one 
may also write the equation 



d Tx d fx T 



^3 {d x d TQ F+Z\ 2m d\T) ?9 3 {d x d T0 F-Z\ 2m d\F) 
®\{Z\2md\T) 



~ {9 f d f l If ^ ) 2 dl log ^3 (Z\ 2m d\T 



(5.30) 



a 2 .f 

e t o 



whose form is close to ()2.30|) but differs by the nontrivial "coefficient" in the square 
brackets. In the limit T — > ioo the theta-function $3 tends to unity, and we obtain the 
dispersionless Toda equation (j2.3()|) . 
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6 Conclusions 



In this paper we have considered the Dirichlet boundary problem in planar multiply- 
connected domains. A planar multiply-connected domain D c is the complex plane with 
several holes. We study how the solution of the Dirichlet problem depends on small 
deformations of boundaries of the holes. 

General properties of such deformations allow us to introduce the quasiclassical tau- 
function associated to the variety of planar multiply-connected domains. By the tau- 
function, we actually mean its logarithm, which only makes sense for the quasiclassical 
or Whitham-type integrable hierarchies. Namely, the key properties are the specific 
"exchange" relations ()3.43|) which follow from the Hadamard variational formula for the 
Green function and the harmonic measure. They have the form of integrability conditions 
and thus ensure the existence of the tau-function. The tau-function corresponds to a par- 
ticular solution of the universal Whitham hierarchy [I] and generalizes the dispersionless 
tau-function which describes deformations of simply-connected domains. 

The algebro-geometric data associated with the multiply-connected geometry include 
a Riemann surface with antiholomorphic involution, the Schottky double of the domain 
D c = C \ D endowed with particular holomorphic coordinates z and z on the two sheets of 
the double, respecting the involution. This Riemann surface has genus g = #{holes} — 
1. The (logarithm of) tau-function, T, describes small deformations of these data as 
functions of an infinite set of independent deformation parameters which are basically 
harmonic moments of the domain. These variables can be equivalently redefined as 
periods of the generating one-form S{z)dz over non-trivial cycles on the double and the 
residues of the one-forms A].(z)S(z)dz, where A^{z) ~ z~ k is some proper global 

basis (jH.lOj) of harmonic functions. 

We have obtained simple expressions for the period matrix, the Abel map and the 
prime form on the Schottky double in terms of the function T . Specifically, all these 
objects are expressed through second order derivatives of the T in its independent vari- 
ables. 

The generalized dispersionless Hirota equations on T for the multiply-connected case 
(equivalent to the Whitham hierarchy) are obtained by incorporating the above men- 
tioned expressions into the Fay identities. As a result, one comes to a series of quite non- 
trivial equations for (second derivatives of) the function F, which have not been written 
before (except for certain relations for the second derivatives of the Seiberg-Witten pre- 
potential [2E])- When the Riemann surface degenerates to the Riemann sphere with two 
marked points, they turn into Hirota equations of the dispersionless Toda hierarchy. 

Algebraic orbits of the universal Whitham hierarchy describe the class of domains 
which can be obtained as conformal images of a "half" of a complex algebraic curve 
with the antiholomorphic involution under conformal maps given by rational functions 
on the curve. In particular, all domains having only a finite amount of non-vanishing 
harmonic moments, are in this class. In this case one can define the curve by a polynomial 
equation written explicitly in [22] ■ This situation is an analog of the (Laurent) polynomial 
conformal maps in the simply-connected case and literally corresponds to multi-support 
solutions of matrix models with polynomial potentials. The definition of the tau-function 
for multiply-connected domains proposed above holds in a broader set-up of general 
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algebraic domains. It does not rely on the finitness of the amount of non-vanishing 
moments. In general any effective way to describe the complex curve associated to 
a multiply-connected domain by a system of polynomial or algebraic equations is not 
known. The curve may be thought of as a spectral curve corresponding to a generic 
finite-gap solution to the Toda lattice hierarchy. 
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Appendix A. Theta functions and Fay identities 

Here we present some definitions and useful formulas from PJ. The Riemann theta 
function 0(W) = 9(W\T) is defined as 

0(W) = ]T e l ™- T - n+2 ^ n W (A.l) 

neZa 

The theta function with a (half-integer) characteristics 8 = (<$i, 8 2 ), where 5 a = T a p8\^ + 
#2,0 and 81,82 £ |Z 5 is 

0^(W) = e in di- T -di+ 27Ti dii w +d2)Q(yv _|_ = ^ e i7r ( n +^i)- T -( n +^i)+ 27ri ( n +^i)-( w +^2) 

n€Z9 

(A.2) 

Under shifts by a period of the lattice, it transforms according to 

e d (w + e a ) = e 2m5 ^e 6 (w) 

9 d (W + T af3 e p ) = e -^, a -^T aa -2,iw ae ^ w ^ 
The prime form E(z,() is defined as 

w( n fl.(W(s)-W(C)) /A 

E(z,Q= , ; — fA4) 



where 9* is any odd theta function, i.e., the theta function with any odd characteristic 
8* (the characteristics is odd if 46* • 8* 2 = odd). The prime form does not depend on the 
particular choice of the odd characteristics. In the denominator, 

30* (W) 



dW a 



w=o 
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is the set of ^-constants. 

The data we use in the main text contain also a distinguished coordinates on a 
Riemann surface: the holomorphic co-ordinates z and z on two different sheets, and we 
do not distinguish, unless it is nessesary between the prime form (|A.4|) and a function 
E(z,Q = E(z, Q{dz) 1 / 2 {dC > ) 1 / 2 "normalized" onto the differentials of distinguished co- 
ordinate. 

Let us now list the Fay identities 0| used in the paper. The basic one is Fay's trisecant 
formula (equation (45) from p. 34 of P) 

0(Wi - W 3 - Z) 9(W 2 - W 4 - Z) E( Zl , z A )E(z 3 , z 2 ) 

+ 0(Wi - W 4 - Z) 9(W 2 - W 3 - Z) E(z u z 3 )E(z 2 , z 4 ) (A.5) 

= e(w 1 + w 2 - w 3 - w 4 - z) e(z) e( Zi , z 2 )e(z 3 , Zi ) 

Here Wj = W(zj). This identity holds for any four points z\,...,Z4 on a Riemann 
surface and any vector Z 6 Jac. In the limit z 3 — > z 4 = oo one gets (formula (38) from 
p. 25 of P) 

^ + /« rfW - Z) 0(Z) S(zx, 2%,) 



- Z) 0(/i 2 rfW - Z) E( Zl , w)E(z 2 , oo) 
dn^)(oo) + f dW a (oo) d z Jo g Yi " 5 



(A.6) 



where 



drt z ^\oo) = d z hg (A. 7) 

rj{z, z 2 J 

is the normalized Abelian differential of the third kind with simple poles at z\ and z 2 
and residues ±1. 

Another relation from [5| we use (see e.g. (29) on p. 20 and (39) on p. 26) is 

9(Wl "^ z ir:; W2+z) -^-^ t (iogs(z)),^ fa) ^M 

U [z 1 ,z 2 ) a f3=1 

(A.8) 

where 

d 2 log#(Z) 

(log0(Z)) >ajB " 



and 

u(zi,z 2 ) = d Zl d Z3 logE(z 1 ,z 2 ) (A. 9) 

is the canonical bi-differential of the second kind with the double pole at z\ = z 2 . 



Appendix B. Degenerate Schottky double 

For an illustrative purpose we would like to adopt some of the above formulas to the 
simplest possible case, which is the Riemann sphere realized as the Schottky double of 
the complement to the disk of radius r. In this case 

dW = — = idip = -— (B.l) 

z z 
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is purely imaginary on the circle and obviously satisfies the condition dW(z)+dW(z) = 0. 
Further, (cf. (ETflHll ) 

W(z) = ( Z dW = j Z dW = log- (B.2) 

and 

w(z) = e w{z) = - (B.3) 

which is nothing but the conformal map of the exterior of the circle | z \ > r onto exterior 
of the unit circle \w\ > 1. Note that on the "lower" sheet of the double 

W{z) = f Z dW = log- (B.4) 

Jr Z 

and instead of (|B.3|) one gets 

w(z) = ~ (B.5) 

which is the conformal map of the exterior of the disk \z\ > r on the lower sheet onto the 
interior of the unit circle \w\ < 1. The prime form on the genus zero Riemann surface is 

(cf. (J23J) 

^/ n Wi — Wo , 

E( Zl ,z 2 )= J_ J- B.6 
V«Wiv dw 2 

where W{ = w(zi). Let us compute E(oo,6b), which is understood in the main text, as 
"normalized" on the values of the local coordinates £00 = 1/2: and z^ = 1/z in the points 
00 and 00 on two sheets of the double. One gets (cf. ([5.9)1 ) 

p2/_ — \ (W(OO) - W(OO)) 2 



Substituting into ()B.7jl the formulas (|B.3|) . ()B.5|1 and 

dz z^ 
dw( 00) = lim — = — lim — c?2; 

Z^OO f z—*oo f 

dw(ob) = r dz^ 

one finally gets 



(B.8) 



z 2 lr 2 1 

£(00, oo) 2 = lim — 1— = - (B.9) 

2— >oo (z r)r t 

and this demonstrates that (j5.14|) indeed turns into (|2.20j) in the limit. 



References 

[1] A. Hurwitz and R. Courant, Vorlesungen iiber allgemeine Funktionentheorie und 
elliptische Funktionen. Herausgegeben und ergdnzt durch einen Abschnitt iiber ge- 
ometrische Funktionentheorie, Springer- Verlag, 1964 (Russian translation, adapted 
by M.A. Evgrafov: Theory of functions, Nauka, Moscow, 1968). 

[2] A. Marshakov, P. Wiegmann and A. Zabrodin, Commun. Math. Phys. 227 (2002) 
131, e-print archive: |hep-th /0109048 



43 



[3] I.M. Krichever, Funct. Anal. Appl. 22 (1989) 200-213 

[4] I.M. Krichever, Commun. Pure. Appl. Math. 47 (1994) 437, e-print archive: 
|hep-th/9205110| 

[5] M. Mineev- Weinstein, P.B. Wiegmann and A. Zabrodin, Phys. Rev. Lett. 84 (2000) 
5106, e-print archive: |nlin. SI/0001007 

[6] P.B. Wiegmann and A. Zabrodin, Commun. Math. Phys. 213 (2000) 523, e-print 
archive: |hep-th /9909147| 

[7] A. Boyarsky, A. Marshakov, O. Ruchayskiy, P. Wiegmann and A. Zabrodin, Phys. 
Lett. B515 (2001) 483-492, e-print archive: |hep-th/0105"260| 

[8] I. Krichever and S. Novikov, Funct. Anal. Appl. 21 (1987) No 2 46-63. 

[9] J.D.Fay, "Theta Functions on Riemann Surfaces", Lect. Notes in Mathematics 352, 
Springer- Verlag, 1973. 

[10] J. Hadamard, Mem. presentes par divers savants a I'Acad. sci., 33 (1908). 

[11] P.J. Davis, The Schwarz function and its applications, The Cams Mathematical 
Monographs, No. 17, The Math. Assotiation of America, Buffalo, N.Y., 1974. 

[12] I. Krichever, 2000, unpublished. 

[13] L. Takhtajan, L ett. Math. Phys. 56 (2001) 181-228, e-print archive: 
|math.QA/0102164| 

[14] I.K. Kostov, I.M. Krichever, M. Mineev- Weinstein, P.B. Wiegmann and A. Zabrodin, 
T-junction for analytic curves, Random matrices and their applications, MSRI pub- 
lications, 40, Cambridge Academic Press, 2001, e-print archive: hep-th/0005259 

[15] E. Hille, Analytic function theory, v.II, Ginn and Company, 1962. 

[16] J. Gibbons and Y. Kodama, Proceedings of NATO ASI "Singular Limits of Disper- 
sive Waves", ed. N. Ercolani, London - New York, Plenum, 1994; 
R. Carroll and Y. Kodama, J. Phys. A: Math. Gen. A28 (1995) 6373. 

[17] K. Takasaki and T. Takebe, Rev. Math. Phys. 7 (1995) 743-808. 

[18] M. Schiffer and D.C. Spencer, Functionals of finite Riemann surfaces, Princeton 
University Press, 1954 

[19] B. Gustafsson, Acta Applicandae Mathematicae 1 (1983) 209-240. 

[20] D. Aharonov and H. Shapiro, J. Anal. Math. 30 (1976) 39-73; 

H. Shapiro, The Schwarz function and its generalization to higher dimensions, Uni- 
versity of Arkansas Lecture Notes in the Mathematical Sciences, Volume 9, W.H. 
Summers, Series Editor, A Wiley-Interscience Publication, John Wiley and Sons, 
1992. 



44 



[21] P. Etingof and A. Varchenko, Why does the boundary of a round drop becomes a 
curve of order four, University Lecture Series, 3, American Mathematical Society, 
Providence, RI, 1992. 

[22] V. Kaza kov and A. Marsh akov, J. Phys. A: Math. Gen. 36 (2003) 3107-3136, e-print 
archive: |hep-th/021l"236| 

[23 

[24 
[25 

[26 

' |hep-th/0306i84| 

[27] A. Marshakov, A. Mironov and A. Morozov, Phys. Lett. B389 (1996) 43-52, e-print 
archive: |hep-th/9607109l 

H. Braden and A. Marshakov, Phys. Lett. B541 (2002) 376-383, e-print archive: 
|hep-th/0205308| 

[28] A. Gorsky, A. Marshako v, A. Mironov an d A. Morozov, Nucl. Phys. B527 (1998) 
690-716, e-print archive: |hep-th/9802007| 



B. De Wit and A. Marshakov, Theor. Math. Phys. 129 (2001) 1504 [Teor. Mat. Fiz. 
129 (2001) 230], e-print archive: |hep-th/0105289l 

F. Gakhov, The boundary value problems, Nauka, Moscow, 1977 (in Russian). 

B.Eynard, Large N expansion of the 2-matrix model, multicut case, e-print archive: 
|math-ph/0307052l 

M.Bertola, Free energy of the two-matrix model/ dToda tau-function, e-print archive: 



45 



